On the size of the cr meson and its nature 
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In this work the nature of the cr or /o(600) resonance is discussed by evaluating its quadratic scalar radius, 
(r 2 )^. This allows one to have a quantitative estimate for the size of this resonance. We obtain that the cr 
resonance is a compact object with (r 1 )^ = (0.19 ± 0.02) - (0.06 ± 0.02) fm 2 . Within our approach, employing 
unitary chiral perturbation theory, the cr is a dynamically generated resonance that stems from the pion-pion 
interactions. Given its small size we conclude that the two pions inside the resonance are merged. A four-quark 
picture is then more appropriate. However, when the pion mass increases, for pion masses somewhat above 
400 MeV, the picture of a two-pion molecule is the appropriate one. The cr is then a spread nn bound state. 
These results are connected with other recent works that support a non standard nature of the cr as well, while 
fulfilling strong QCD constraints, as well as with lattice QCD. 

We offer a detailed study of the low-energy S -wave mx scattering data from where we extract our values for 
the threshold parameters of S -wave mx phase shifts, the 0(p 4 ) chiral perturbation theory low energy constants as 
well as the cr pole position. From the comparison with other accurate determinations in the literature we obtain 
the average values for the isospin S -wave nn threshold parameters, a° = 0.220±0.003, b% = 0.279±0.003 M~ 2 , 
and for the real and imaginary parts of the cr pole position in ^fs, 458 ± 14 - j 261 ± 17 MeV. The quark mass 
dependence of the size of the cr, its mass and width are considered too. The latter agree accurately with a 
previous lattice QCD calculation. The fact that the mass of this resonance tends to follow the threshold of two 
pions is a clear indication that the cr is a dynamically generated meson-meson resonance. 

PACS numbers: ll.10.St, 12.39.Fe, 11.30.Rd, 14.40.-n,14.40.Be, 14.65.Bt, 13.75.Lb 

Keywords: Chiral Symmetry; Non-perturbative methods; Sigma meson; Scalar radius; Hadron size; Form factor; Lattice QCD. 



I. INTRODUCTION 

The lightest resonance in QCD with the quantum numbers 
of the vacuum, J PC = ++ , is the cr or /o(600) resonance [1]. 
Its connection with chiral symmetry has been stressed since 
the sixties in the linear sigma model [2], while its tight relation 
with the non-linear sigma model was realized in the nineties. 
In this respect there have been several papers that clearly con- 
nect this resonance with chiral dynamics of the two-pion sys- 
tem. One has first to mention the works of Truong and collab- 
orators [3-6] who first emphasized the important role played 
by the null isospin (I) S -wave nn final state interactions in 
several processes giving rise to a strong numerical impact on 
the estimations based on current algebra technique or chiral 
perturbation theory (ChPT) [7-11]. A notoriously improved 
comparison with experiment was then obtained, e.g. for K( A 
decays [3], 77 — » 37r [4], scalar and pion vector form factors 
[5] and nn scattering [6]. These works stress the role of the 
right-hand or unitarity cut and make use of a method to resum 
unitarity based on the expansion of the inverse of a form factor 
or scattering amplitude. This is the so called Inverse Ampli- 
tude Method (IAM), that in the end is analogous to a Padde 
method of resummation. Within this technique the cr pole was 
first obtained in Ref. [12], together with the K* and p reso- 
nances in the P-waves. However, due to the lack of coupled 
channels, no further light scalar resonances were generated, 
in particular the /o(980) and a (980). 

Independently, the cr resonance pole was also obtained si- 
multaneously in Ref. [13], together with the / = /o(980) 



* albaladejo@um.es 
' oller@um.es 



and ao(980) resonances. The associated amplitudes were de- 
termined by solving the Bethe-Salpeter equation taking as po- 
tential the lowest order ChPT Lagrangian. Only one free pa- 
rameter (a natural sized cut-off) was involved. Later on, when 
the IAM was extended to coupled channels [14], it was pos- 
sible to obtain in Refs. [14-16] the cr, / (980), a (980) and 
k resonances altogether, that is, the whole nonet of the light- 
est scalar resonances [17-22], together with the nonet of the 
lightest vector resonances. 

The approach of Ref. [13], based on solving a Bethe- 
Salpeter equation, was put on more general grounds in 
Ref. [19] by applying the N/D method [23]. In this way, it 
is possible to include higher orders in the chiral counting 
[24, 25] as well as explicit resonant fields [26], if required. 
Later works based on this scheme are Refs. [27-30]. With this 
approach [19] one builds a unitarized meson-meson scattering 
amplitude by solving the N/D equation in an algebraic way so 
that an approximate solution is obtained by treating perturba- 
tively the crossed cuts. As a result, the ChPT expansion is re- 
produced order by order, while the unitarity cut is resummed 
[24]. In this respect, one should stress that the crossed cuts 
can be treated perturbatively for the isoscalar nn S -wave. Its 
size was estimated to be smaller than 10% in Ref. [19] along 
the physical region for energies up to around 1 GeV. Indeed, 
different approaches with various degrees of sophistication 
provide very similar values for the cr pole resonance param- 
eters, mass and width. Either by employing just the leading 
order (LO) ChPT [13] (without left-hand cut at all), next- 
to-leading order (NLO) [12] or next-to-next-to-leading order 
(N 2 LO) [31]. In these two later references the left-hand cut is 
included as calculated by ChPT at one and two-loop orders, 
respectively. The fact that the results are very similar clearly 
indicates that the left-hand cut is indeed a perturbation. The 
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cr pole positions in yfs, with s the total center of mass energy 
squared, obtained in these works are: yfs^ = 468 - ; 194 MeV 
[13], 440 - i 245 MeV [12] and 445 - i 235 MeV [31]. In 
the following we identify the mass and half width of the cr 
resonance from the pole position as M a — i rV/2 = y[s^. 

More recently, Ref. [32], based on the solution of the Roy 
equations [33] and ChPT at two-loops [34, 35], obtained 
the value 445^ 6 - i 272+^ MeV. The Roy equations im- 
plement crossing symmetry exactly, while the previous ref- 
erences [12, 13, 19, 31] do it perturbatively. The fact that all 
these pole positions for the cr lie rather close to each other 
(particularly one can say that convergence is reached very ac- 
curately for the real part) is another indication for the correct- 
ness of treating crossed-channel dynamics perturbatively, as 
done in the framework of Refs. [19, 36] (see also [37-39]). 
Indeed, to our mind, both schemes are complementary be- 
cause the Roy equations need for their implementation of the 
knowledge of large amount of data in several partial waves 
up to high energies, which is affected by systematics errors 
in experiments (many of them old ones) and also in theory 
{e.g. high energy extrapolations), not always well under con- 
trol. Let us also mention that all these analyses neglect alto- 
gether the inelasticity due to the An channel in nn S -waves so 
that, up to the opening of the KK threshold at around 1 GeV, 
no inelasticity is assumed. The An channel was approached 
in Ref. [28] as crcr and pp states (with their couplings to all 
the channels predicted from chiral dynamics) and found the cr 
pole at 456 ± 6 - i 241 + 7 MeV. 1 This pole position is quite 
close to those in the previous references and compatible with 
the result 484 ± 17 - / 255 ± 10 MeV from Ref. [40]. Thus, 
since the pole positions of Refs. [12, 13, 28, 31, 32, 40] lie so 
close to each other we could conclude that our present knowl- 
edge on the pole position of the cr resonance is quite precise 
and, furthermore, we understand the underlying physics at the 
hadronic level. 

Between earlier approaches to the previous discussed re- 
sults based on ChPT concerning the lightest scalars, we have 
Refs. [17, 41] within the MIT bag model that already in the 
late seventies predicted a complete nonet of four-quark ++ 
resonances (comprising the cr, /o(980), «o(980) and k), with 
M a = 660 + 100 MeV and IV = 640 ± 140 MeV. The four- 
quark nature of the lightest scalars is also favored in Refs. [42- 
45] attending to scattering and production data, including two- 
photon fusion, J/*¥ and cf> decays, and in Refs. [46, 47]. The 
important role played by two-meson unitarity for understand- 
ing the scalar sector for yfs < 1 GeV was also stressed in 
Ref. [18] (a similar approach was later followed in Ref. [48]), 
employing a unitarized chiral quark model, and in Ref. [49], 
within the Jiilich meson-exchange models. Considerations 
based on increasing the QCD number of colors, Nc, were ex- 
ploited in Refs. [19, 47, 50-52], showing that the cr resonance 
has a non-standard Nc dependence. This can be done more 
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safely for Nc ^ 3, not too large, while statements for Nc » 3 
depend much more on fine details of the approach [39, 52-57]. 
QCD sum rules were also applied for the study of the lightest 
scalar meson, e.g. in Refs. [58-62]. It is argued too that the 
cr resonance is the chiral partner of the pion [63-66] and the 
way in which the cr pole evolves when approaching the chiral 
symmetry restoration limit is different according to the nature 
of this resonance [67]. 

From an experimental point of view new interest is trig- 
gered on the cr resonance from recent high-statistics results, 
e.g. J/^V — > conn where a conspicuous peak is seen [68]. In- 
deed, this decay mode was the first clear experimental sig- 
nal of a cr resonance [69, 70]. Another marked peak around 
the cr energy region is also observed in several heavy me- 
son decays. E.g. it was observed with high statistical signif- 
icance in D — » n + n~n + [71]. Both types of decays present a 
strong peak in the cr mass energy region because the absence 
of the Adler zero in the pion scalar form factor, as explained 
in Refs. [72, 73]. 2 Another field of increasing activity, both 
experimental [76, 77] and theoretical, concerns the fusion of 
two photons into a pair of pions and from there to extract the 
width of the cr to yy [78-82]. This is also expected to shed 
light on the nature of the cr meson [78]. 

The relative strength of the cr coupling to KK compared to 
nn is also taken as an important property in order to disen- 
tangle between different models for the nature of the cr me- 
son {qq, four-quarks, glueball or ^-molecule), as stressed in 
Ref. [59]. This reference points out that the not so much sup- 
pressed coupling of the cr to K + K~ (gaK+K-), as compared 
with that to n + n- (g a7l+7l -), l&r/r aH/I W*- 1 = 0.37 + 0.06 
[59], is a key ingredient to advocate for a gluonium nature 
of the cr meson. According to Ref. [59], a simple qq inter- 
pretation of the cr fails to explain the large width of the cr 
while a four-quark scenario has difficulties to explain its large 
coupling to K + K~. It is then worth emphasizing that the T- 
matrices obtained in Refs. [13, 19, 28] also predict a ratio 
for the cr couplings to K + K~ and n + n~ in perfect agreement 
with the value above of Refs. [58, 59, 83]. Explicitly, we have 
\ga-K*K-\/\gan*7r\ = 0.36 ± 0.04 from the average collected in 
Ref. [20]. However, in our case this stems from the dynamical 
generation of the cr resonance from the Goldstone boson dy- 
namics associated to the strong scalar isoscalar nn interaction. 
We also stress that this approach has been confronted with a 
large amount of data from different reactions, both scattering 
and production experiments, in most of the reactions already 
quoted in this introduction. 

One of the aims of this work is to show that the often iden- 
tification of dynamically generated resonances from the inter- 
actions of two mesons (pions in our case) as meson-meson 
molecules is misleading. As we show here, depending on the 
meson mass, one can have situations where the size of a dy- 
namically generated meson-meson resonance is certainly too 
small to be qualified as a two-meson molecule. Indeed, its size 
could be as small as that of one of the mesons involved in their 
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formation. The fact that the <x is such a tight compact object 
clearly hints that the two pions pack so much that it is not 
meaningful anymore to keep their identities separately. At this 
stage, a four quark compact resonance seems a more appropri- 
ate picture. This is also supported by the Nc evolution of the 
cr-pole trajectory which is clearly at odds with the expecta- 
tions for a purely qq or glueball resonance, but in the lines 
of what it is expected for a meson-meson or four quark reso- 
nance [39, 52-57, 84]. However, by increasing the pion mass 
the cr resonance pole tends to follow the two pion threshold, 
and when it is close to the latter its size increases, becoming a 
spread object. This is a clear indication for the molecular char- 
acter of the cr for large enough pion masses, M„ > 400 MeV. 
In addition, let us also emphasize that our work is the first cal- 
culation of the size of the cr resonance. This is a novel way to 
study its nature in the literature. 

The rest of the paper is organized as follows. In Sec. II we 
give a short introduction to the S U (2) ChPT Lagrangians both 
at LO and NLO that are used in the rest of the paper. Next we 
dedicate Sec. Ill to evaluate nn scattering at one-loop order in 
Unitary ChPT. A wide set of data is fitted, including some re- 
cent lattice QCD determinations as a function of M„. We pay 
special attention to the threshold parameters and the cr pole 
position. For these quantities we also compare with previous 
phenomenological determinations and the lattice QCD results 
on the dependence of the cr pole mass as a function of the pion 
mass. We dedicate Sec. IV to the calculation of the scalar form 
factor of the cr resonance. First pion scattering in the presence 
of a scalar source is discussed. The scalar form factor of the cr 
is calculated from the double cr pole present in the amplitude 
for the previous process, once nn initial and final state interac- 
tions are taken into account. Then, we determine the quadratic 
scalar radius of the cr and then have some information on the 
size of this resonance. We stress that this radius is pretty small, 
around 0.5 fm indicating that the cr is a compact object. We 
also discuss the relation between the value of the cr scalar form 
factor at the origin and the dependence of the cr pole with the 
pion mass, related by the Feynman-Hellmann theorem. Both 
issues, the quadratic scalar radius and the Feynman-Hellman 
theorem, are addressed in Sec. V. After concluding in Sec. VI, 
we dedicate Appendix A to the loop functions used through- 
out the amplitudes calculated, which are in turn given in Ap- 
pendix B for pion scattering in the presence of a scalar source. 



II. S U(2) CHIRAL LAGRANGIANS 

We follow the standard ChPT counting and the processes 
under consideration, the scattering of pions with and without 
the presence of a c-number external scalar source, are calcu- 
lated both at LO and NLO. The chiral power counting of a 
connected diagram, p D (where p is a generic small momen- 
tum compared to AchPT - 1 GeV), obeys the equation [7, 10] 



D = 2 + Y N d (d -2) + 2L 



(1) 



loops. Each derivative increases the counting by one unit and 
the lightest quark masses add two units to D. The LO calcu- 
lation has D = 2 with no loops (L - 0) and involves only 
d — 2 vertices. For the NLO, D — 4, and one has diagrams 
with L — 1 that involve only d — 2 vertices. There are also 
diagrams with L — with only one d = 4 vertex, with the rest 
of vertices having d = 2. 

Up to NLO, 0(p 4 ), one has to consider the S 1/(2) chiral 
Lagrangians at 0(p 2 ), £a, and 0(p 4 ), JU, that we take from 
Ref. [7]: 

£2 = y(V„tf T V*tf + 2&r T l/)) , (2) 

£ 4 = Zj (wt/ T v M t/) 2 + 1 2 (vc/ T v v t/) (v^VyU) 

+ h( X J uf + U(V X T V,U) + --- (3) 

where F is the pion weak-decay constant in the chiral limit, 
the terms proportional to the Z, are the NLO chiral countert- 
erms and the ellipsis indicate terms not shown because are not 
needed here. The pion fields are included through the 0(4) 
real vector field U(x) of unit length, U T U = 1, as: 



U T = (u°(x), 0(x)) , 



U\x> 



1 



1 



U(x) 2 = 1 - -U L - -ET- 



(4) 



We also use the relation between the charged and Cartesian 
pion fields given by 

1 — • 2 

7T 4- ITT „ t 

, n°=n\ (5) 



V2 

The explicit chiral symmetry breaking due to the finite u 
and d quark masses enters through the vector-field x T ( x ) = 
2B(m + s(x), p'(x)). Here, 2Bm = M 2 is the pion mass at lead- 
ing chiral order 3 and m is the algebraic mean of the u and 
d quark masses (we consider exact isospin symmetry). The 
fields s(x) and p'(x) refer to the scalar and pseudoscalar c- 
number external sources, in order. The covariant derivative 
reduces in the problem that we are studying to the standard 
derivative, — > d^, since we do not consider here external 
vector nor axial-vector currents. Finally, the parameter B is re- 
lated to the value of the quark condensate in the chiral limit 
(q'qi) = -S ij F 2 B [7]. 

In the following we employ the finite and scale independent 
constants Zj defined by 

k = i r i+yi- R 



R 



l 32n 2 ' 

M 2 



y, 



,n-4 



(io g 4^+r'(i) + i) 



(6) 



In this equation, d is the chiral dimension of a vertex, Nd the 
number of vertices with dimension d and L is the number of 



3 In our notation M represents the pion mass at LO, i.e., the parameter that 
appears directly from the Lagrangian, while M„ refers to the physical pion 
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FIG. 1 . Diagrams for the one-loop calculation of the pion self-energy. 
Full circles represent 0{p~) vertices, while the empty ones corre- 
spond to the 0(p 4 ) vertices. 



so that li is, up to a numerical factor, the renormalized cou- 
pling constant l r . at the scale p = M M p . In the chiral limit 
the li are not defined as they are then divergent quantities. The 
needed y,- coefficients are [7]: 



12 1 

7i = j , 72 = 3 , 73 = . 74 = 2 



(7) 



The infinite quantity R is cancelled with the infinities that orig- 
inate from loops, see Appendix A. 

The calculation of the pion self-energy, -iL(p 2 ), is neces- 
sary in order to take into account the renormalization of the 
wave function of the initial and final pions. One has: 



Up 2 ) = 



3M 2 A (M 2 ) 2M 4 h p 2 A Q (M 2 ) 



2F 2 



F 2 



F 2 



(8) 



Notice that Up 1 ) is linear in its argument. The one-point func- 
tion Ao(M 2 ) is given in Eq. (A2), Appendix A, together with 
the different n-point loop function needed in this work. We 
can write the self-energy Eq. (8) as: 



Z(p 2 ) = Z(M 2 ) + Z(Ml)(p 2 -M 2 ). 



(9) 



The Dyson resummation gives for the renormalized propaga- 
tor, A R (p 2 ), 



i&Rip 2 ) 



iZ 



(p 2 -Ml){\-I.'{Ml)) p 2 -M 2 K 



(10) 



where 



Z 1 + 5Z = 1 



A (M 2 ) 



F 2 



+ 0{M A X ), 



Mt = Af 1 - 



M 1 



32n 2 F 2 



h\+0{M b n ). 



(11) 



Then, in order to take into account the renormalization of the 
pion wave function in our diagrams (both for nn scattering and 
for the nns — > 7T7t process), with four external legs, we have 
to multiply by a factor (Z 1/2 ) 4 = Z 2 = 1 + 26Z + 0(M 4 ). In 
the following, we should keep in mind that the pion propaga- 
tors employed are /A«(p 2 ), Eq. (10), in terms of the physical 
pion mass. This will make simpler the calculation of some di- 
agrams for the process nns — > nn. Let us also mention that the 
amplitudes calculated are given in terms of the physical mass 
and weak decay constant of the pion. The latter is given by 
[7]: 



F„ = F 1 



M 2 



\6n 2 F 2 



h\ + 0{Ml) 



(12) 




(c) 




FIG. 2. Feynman diagrams for mr scattering up to NLO. Full cir- 
cles represent 0(p 2 ) vertices, while the empty ones correspond to the 
0(p 4 ) vertices. 



HI. nn SCATTERING AND THE a MESON 
A. The nn —> nn amplitude 

The chiral Lagrangians exposed in Sec. II comprise four 
low energy constants (LECs), at 0(p 4 ). Additionally, our 
resummation procedure, explained below, includes a subtrac- 
tion constant through the two-meson unitarity one-loop func- 
tion. Before considering the nns — > nn amplitude, we must 
fix these free parameters. This is accomplished by comparing 
our results for the scalar nn — > nn phase shifts with / = 0, 2 
with experiment, and also other observables with lattice QCD 
determinations. 

We denote by Xnis, f) the / = nn scattering amplitudes 
calculated from Fig. 2 in ChPT at 0(p"), with n — 2 or 4. 
Their projection in S -wave are indicated by f„( s). Diagram a) 
is the LO contribution, while the rest of diagrams are the NLO 
ones. The last two diagrams, namely, e) and f) contribute to 
the wave-function renormalization of the pion external legs. 
We introduce the usual Mandelstam variables s, t and u. The 
variable s corresponds to the total energy squared of the two 
pions in their center of mass frame (CM), while the other two 
are defined as: 

t = -2p 2 (l -cos0) 
u = -2p 2 (l + cosfl) 



s + t + u = AM 2 

P 2 ^-M 2 



(13) 



Here, p is the three-momentum squared of the pions in their 
CM and 9 is the scattering angle in the same reference frame. 
The amplitudes %„{s) are then given by, 



£.(*) 



i r 1 



(14) 



In the previous equation an extra factor of 1 /2 has been in- 
cluded, in correspondence with the so called unitarity normal- 
ization [13]. The / = nn state is symmetric under the ex- 
change of the two pions so that the unitarity normalization 
avoids having to take into account the presence of the factor 
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1/2 whenever it appears as an intermediate state. In this way, 
the same formulas as for distinguishable particles can be em- 
ployed. In the following of the paper we employ the unitarity 
normalization in all the isoscalar nn matrix elements unless 
the opposite is stated. 

Let us indicate by T(s) the scalar-isoscalar unitarized nn 
partial-wave amplitude. Following the unitarization method of 
Refs. [19, 24], the right-hand cut or unitarity cut is resummed 
by the master formula: 



T{s) = 



V(s) 



1 + V(s)G(s) 



(15) 



This formula is deduced by solving algebraically the N/D 
method [19, 23], treating perturbatively the crossed cuts, 
whereas the unitarity cut is resummed exactly. Here, G(s) is 
the scalar two-point function, 



1 / M 2 „ a-(s) - 1 

G(i) = m\ a + lo § ~r - °"( J ) lo g "TTT7 

\6n z \ p 1 cr(s) + 1 



(16) 



with chiral order p . In the previous equation <r(s) = 
-Jl - 4M 2 /(s + ie). The interaction kernel V(s) has a chiral 
expansion, V{s) = V2(s) + Vi(s) + ■ ■ ■ , with the chiral or- 
der determined by the subscript. The different chiral orders 
of V(s) are calculated by matching T(s) with its perturbative 
expansion calculated in ChPT In this way up to 0(p 4 ), 



T(s) 



V(s) 



1 + V(s)G(s) 
= fr(s) + fr(s) + ■■■ 
= V 2 (s) + V 4 (s) - V*(s)G(s) + . . 



(17) 



where the ellipsis indicate 0(p 6 ) and higher orders in the ex- 
pansion. It results then: 



V 2 (s) = , 
V A (s)={ A (s)+t z (s) 2 G(s). 



(18) 



The finite piece of the unitarity term in Fig. 2 (that is, the 
term of £t(s) that contains the unitarity cut and is proportional 
to the unitarity two-point one-loop function) is given by: 



-£(s)B (s) 



(19) 



Here, Bq(s) is the two-meson loop in dimensional regulariza- 
tion, without the R + log(M 2 / p 2 ) piece (that cancels out with 
the other infinite and scale dependent terms, see Eqs. (A5) 
and (A6) in Appendix A). In this way, the kernel V{s) = 
V 2 (s) + V4(s) has no unitarity cut because: 

+ fi(s)G(s) = -f 2 (s)(Bo(s) - G(s)) , (20) 

and the cut cancels in the r.h.s. of the previous equation. The 
full unitarity cut arises from the denominator 1 + V(s)G(s) in 
Eq. (15). 

In this Section we have dealt with the 1 = unitarized am- 
plitudes but, needless to say, the same formalism applies to the 
1 = 2 ones, by just changing the kernel V(s). We additionally 
note here that the same subtraction constant is used for both 
channels, as required by isospin symmetry [85]. 



TABLE I. Summary of our LO and NLO fits. In the last column the 

x 2 



' 2 per degree of freedom is given. 



Fit 


a 




x 1 

rl.o.f. 


LO 


-1.36 ± 


0.12 


1.6 


NLO 


-1.2 ± 


0.4 0.1 


! ± 0.9 4.6 ±0.4 2 + 4 3.9 ± 0.5 0.7 



B. Fits and the cr meson 

At LO, there is just one free parameter corresponding to 
the subtraction constant in G(s). At NLO, there are, in addi- 
tion, four LECs, /,, i = 1,2,3,4. For 1 = 0, the phase shifts 
that we fit contain the very precise data of K e \ decays be- 
low = 400 MeV [86-90]. These data are corrected for 
isospin breaking effects, as explained in Ref. [91]. Above that 
energy, the data of Ref. [92] and the average of different ex- 
periments [93-98], as used e.g. in Ref. [19], are taken into 
account. For 1 = 2, the data come from Refs. [99, 100]. The 
fits extended to a maximum energy y/s m:lx = 0.8 GeV at LO, 
both for 1 = and 1 = 2, whereas at NLO we extend this 
range up to ys max = 1 GeV for 1 = 2. This is not done for 
1 = because of the related presence of the KK threshold 
and the /o(980) resonance. The phase shifts are denoted by 
6q, with / = 0, 2. For our NLO fits we also fit recent lattice 
QCD results as functions of the pion mass for F„ [101, 102] 
and the isotensor scalar scattering length, a? [102, 103]. 4 The 
dependence of F„ with the pion mass is calculated at NLO in 
ChPT, Eq. (12). The scattering length is defined through 
the threshold expansion in powers of p 2 of our full results: 



fc uP 2 + 0(lpl 4 ), 



(21) 



that we extrapolate in terms of the pion mass squared. 

The resulting values for the fitted parameters are given in 
Table I. At LO the subtraction constant for the G(s) function 
is a = -1.36 ±0.1 2. Four LECs appear additionally to the sub- 
traction constant as free parameters at NLO. In order to avoid 
large correlation among them, the subtraction constant at NLO 
is constrained to remain near its value at LO. This is done by 
adding a new term to the y 1 taking into account the differ- 
ence between the values of a at LO and NLO, but enlarging 
its error at LO from 0.12 to 0.2, so that its contribution to the 
resulting^- 2 is tiny but enough to remove the large correlations 
that would appear otherwise among the LECs and the subtrac- 
tion constant. The parameters of both fits (LO and NLO) are 
shown in Table I, and the corresponding phase shifts are plot- 
ted in Fig. 3 with their respective errors. The left panel is for 
1 = and the right one for 1 = 2. The (red) dashed lines arise 
from our fit at LO (V(s) = V 2 {s)), whereas the (blue) solid 
ones show the NLO fit (V(s) = V 2 (s) + V^s)). In the inset of 



4 We consider the spread of these lattice QCD results as a source of system- 
atic error for our fits. The final errors included in the fit are depicted by the 
dashed error bars in Fig. 4. 
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FIG. 3. Comparison of our scalar nn phase shifts to experimental data for / = (left panel) and 1 = 2 (right panel). The (red) dashed line 
shows our fit for the LO case (V(s) = Vi(,s)), whereas the (blue) solid one shows the NLO fit (V(s) = Vifa) + V4(*)). The bands represent our 
uncertainties. The inset in the left panel shows in more detail the low energy K e 4 decays data. The data for / = are from the K e 4 decay data 
of Refs. [86-90] (with isospin breaking effects taken into account as in [91]) and other data from Refs. [92-98]. For 1 = 2 the phase shifts are 
fromRefs. [99, 100]. 
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FIG. 4. Dependence of F n (left panel) and (right panel) with M n as compared with lattice QCD data. The (blue) solid line is given by our 
NLO fit, whereas the band represents our estimated error. The data are from Refs. [101-103]. For we also show the data of Ref. [104], 
although we do not include them in our fits. 



the upper panel the agreement of our results with the lower en- 
ergy data from K e 4 decay can be appreciated. We must stress 
that the difference between LO and NLO manifests mostly in 
the 1-2 channel phase shifts, as can be seen in Fig. 3. In this 
channel, the left-hand cut is more important, but our ampli- 
tudes only incorporates the latter in a perturbative way, so that 
at NLO it is well reproduced, but it is absent at LO. In Fig. 4 
our results for F„ (left panel) and (right panel) are shown, 
and compared with the aforementioned lattice QCD results. 

In Table II we collect some phenomenological [7, 34, 105- 
1 10] and lattice QCD [111-1 16] determinations of the LECs. 
For the latter the last values of each collaboration are taken, 
and, in addition, the direct SU(2) fit results are selected 
if values for SU(2) and SU(3) fits are offered. We have 
also included the range obtained for 1$ from the data of the 
NA48/2 Collaboration [90]. These determinations are com- 
pared graphically in Fig. 5, where for every LEC the different 
results are compatible within errors. The lattice QCD results 
concerning l\n are scarce. The JLQCD and TWQCD Collab- 
orations [114] recently reported l\ -l 2 = -2.9 + 0.9 + 1.3, 
whereas, from our fit, we obtain l\ — h = -3.8 + 1.3. For the 



phenomenological determinations in Table II, since l\ % agree 
well between each other, also the aforementioned difference 
between these LECs does. We finally note that from our fit 
we obtain at NLO ChPT that F = 86.8 + 0.8 MeV, so that 
F„/F = 1 .065 + 0.010, compatible with the estimate of lattice 
QCD results given in Ref. [117], FJF = 1.073 ± 0.015. 

Our function G(s) stems from the calculation of a once- 
subtracted dispersion relation (see e.g. Ref. [19]). If, instead, it 
is calculated approximately by employing a three-momentum 
cut-off A, one has the following relation between the subtrac- 
tion constant and A [15, 24]: 

eu 2 (M 2 \ 
«(/,) = -l+log^+O^fj . (22) 

Our values for the fitted subtraction constant gives a cut-off 
A 750 MeV M p , which is quite a natural value. We will 
make use of these considerations based on Eq. (22) later on, 
when dealing with the M„ dependence of the cr pole position. 

The cr pole appears in the second or unphysical Riemann 
sheet of the amplitude. This sheet is reached by changing the 
function G{s) in the following manner [13]. For s real and 
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TABLE II. Comparison of different phenomenological and lattice QCD determinations of the LECs ( = 1,2,3,4. Together with every 
reference, for an easier comparison the initials of the authors or those of the collaboration are given. 



Ref. 




h 




h 


h 




k 


[7] GL 


-2.3 ± 


3.7 


6.0 ± 


1.3 


2.9 ± 2.4 




4.6 ±0.9 


[34] CGL 


-0.4 + 


0.6 


4.31 


t0.ll 


- 




4.4 ± 0.2 


[105] ABT 


0.4 ± 


2.4 


4.9 ± 


1.0 


2 5 +1 ' 9 




4.20 ±0.18 


[106] PP 


-0.3 ± 


1.1 


4.5 ± 


0.5 


- 




- 


[107] GKMS 


0.37 


±0.95 ± 1.71 


4.17 


±0.19 + 0.43 


- 




- 


[108] BCT 




- 




- 


- 




4.4 ±0.3 


[109] OR 




- 




- 






4.5 ± 0.3 


[110] DFGS 




- 




- 


-15 ±16 




4.2 ± 1.0 


[90] NA48/2 










2.6 ±3.2 






[111] RBC/UKQCD 










2.57 + 0.18 




3.83 ±0.9 


[112] PACS-CS 










3.14 + 0.23 




4.04 ±0.19 


[113] ETM 










3.70 ± 0.07 


±0.26 


4.67 ±0.03 ±0.1 


[114, 115] JLQCD/TWQCD 










3.38 + 0.40 


± 0.24+JJ- 31 


4.09 ±0.50 ±0.52 


[116] MILC 










2.85 ±0.81^ 


3.98 ± 0.32+0^ 


This work 


0.8 ± 


0.9 


4.6 ± 


0.4 


2 + 4 




3.9 ±0.5 



This work 

GKMS 

PP 

ABT 

CGL 

GL 



3.5 



This work 
NA48/2 

JLQCD/TWQCD 
ETM 
PACS-CS 
MILC 

RBC/UKQCD 

ABT 

GL 



This work 

GKMS 

PP 

ABT 

CGL 

GL 
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OR 

DFGS 

BCT 

ABT 

CGL 

GL 



3.5 



4.5 
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FIG. 5. Comparison of the different lattice QCD and phenomenological determinations of the LECs collected in Table II. The (green) diamonds 
are lattice QCD determinations, and (red) circles are the phenomenological ones. The range obtained for /3 by the NA48/2 Collaboration is 
represented by a (blue) triangle. The (black) squares are our results. For an easier comparison, we have included a shaded area that represents 
our results (except for [3). 



above threshold we have 

GnO + ie) = d(s + ie) - AG(s) , (23) 

where the subscript denotes the physical (I) or the unphysi- 
cal (II) Riemann sheet. In the previous equation, AG(s) is the 



discontinuity along the unitarity cut, 

AG(s) = GiO + ie) - Gi{s - ie) = -i^= , (24) 

87r ys 

with p(s) = ^Jlp = yjs/4 - M%, the CM pion three- 
momentum, such that Im p(s) > 0. In order to explore the 
unphysical Riemann sheet, one then makes the analytical ex- 



TABLE III. cr pole position and threshold parameters for the 
isoscalar scalar partial-wave amplitude. 



Fit 


Vv (MeV) 


a 
"o 


b°M 2 


LO 465 


± 2-f231± 7 


0.209 ± 0.002 


0.278 ± 0.005 


NLO 440 


+ 10-; 238 ± 10 


0.219 ± 0.005 


0.281 ±0.006 



0.29 



0.285 



* 0.28 



0.275 





0.225 



445 450 455 
Rey^ (MeV) 



470 



FIG. 6. Montecarlo-like error analysis for the cr mass (M^ = Re -^v) 
and half-width (T^/2 = -Im ys~x) and threshold parameters aJJ and 
bg. The (blue) error ellipses correspond to the NLO fit while the sin- 
gle (red) point with errors is for the LO result. 



trapolation in the cut complex s plane of Eq. (23). 

In the second sheet the cr resonance is a pole in the 1 = 
S -wave 7T7T amplitude, 



Tn(s - So-) = — 



(25) 



being go- the coupling to the nn channel and the ellipsis in- 
dicate the rest of terms in the Laurent series around v (with 
Imv < 0). The pole position v is given in Table III, together 
with the resulting values for the threshold parameters of the 
scalar-isoscalar partial wave. The cr pole position is used to 
define its mass and width, M a - zIV/2 = yfs^. 

The error analysis for any quantities calculated here (e.g. 
the fitted values for the LECs, cr pole position, etc) is per- 
formed by randomly varying our parameters around their fit- 
ted values and accepting those values for the parameters which 
have a v 2 < v 2 . + Ay 2 . Here v 2 ■ is the best value for the v 2 . 

A A mm A A- mm A 

For the LO case, since there is just one free parameter, we give 
our two-sigma confidence interval (otherwise the errors would 



TABLE IV. Values of r^/Z, a° and b° extracted from the litera- 
ture. The value of Ref. [90] corresponds to the latest experiment on 
K e 4 decays (with the errors added in quadrature for an easier com- 



parison). 


Ket. 


M a (MeV) T CT /2 (MeV) 


< 




[118] 


470 ± 50 


285 ± 25 








[32] 


441+' 6 


272+^3 








[40] 


484 ± 17 






±0.013 


0.285 ±0.012 


[28] 


456 ± 12 


241 ± 14 








[119] 


463 ± 6*1) 


254 ± 6 + _H 


0.218 


±0.014 0.276 ±0.013 


[59] 


452 ± 12 


260 ± 15 








[120] 


457+]^ 


279+] 1 








[90] 






0.222 


±0.014 




[35] 






0.220 


± 0.005 


0.276 ± 0.006 


This work 440 ± 10 


238 ± 10 


0.219 


± 0.005 


0.281 ±0.006 


Average 


453 ± 5 


258 ± 5 


0.220 


± 0.003 


0.279 ± 0.003 


Mean 


458 ± 14 


261 ± 17 


0.223 


± 0.007 


0.280 ± 0.004 



be too small), given by A^ 2 = 4. At NLO the one-sigma con- 
fidence interval corresponds to Ay 2 = 5.9. The resulting error 
ellipses are shown in Fig. 6 for the threshold parameters, up- 
per panel, and for the cr mass and width, lower panel. Notice 
that since there is only one free parameter at LO then a curve 
results instead of an error ellipse as in NLO. This is why at 
LO we have just shown the resulting value with its errors. 



C. The cr meson. Comparison with other determinations 

We compare now our results for the cr mass and width as 
well as for the threshold parameters with other determinations 
from Refs. [28, 32, 40, 59, 118-120]. References [32, 40] 
are recent sophisticated determinations of the pion pole po- 
sition claiming to be very precise. In Ref. [28], based on chi- 
ral Lagrangians and the implementation of the N/D method, 
a detailed study of meson-meson scattering in the scalar sec- 
tor up to around -\fs — 2 GeV was performed. All the rele- 
vant channels were taken into account, even the 47r channel 
through the crcr and pp channels whose interactions kernels 
were predicted making use of chiral symmetry and vector me- 
son dominance. A good description of the data considered was 
achieved, which allowed a full description of the resonances 
experimentally seen up to that energy. 5 

The relevant quantities contained in those references are 
collected in Table IV, and compared in Figs. 7 and 8 with 
our LO and NLO determinations. If all these determinations 
can be considered as different measures of the same physi- 
cal quantity, then they should be compatible. A good check 
of their mutual compatibility is to determine whether they are 



5 In Table IV we double the errors of our previous determination [28], so as 
they have a similar size as those from other calculations. In this way the 
weighted average is not so much biassed from just one determination. 
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FIG. 7. In this figure we show the values for the threshold parameters ajj and fcjj from different papers in the literature, as indicated in the plots. 
In the first two panels, from left to right, the (dark gray) inner strip corresponds to the interval covered by the weighted average whereas the 
(light gray) outer strip is for the mean value, both given in Table IV. In the last panel, the rectangles correspond to the aforementioned intervals 
in the a^-b° Q M^ plane. The references are: A [ 1 19], B [40], C [35] and D refers to the NLO determination of this work. 



420 440 460 480 
M a (MeV) 



500 220 240 260 280 
T a /2 (MeV) 



B 

C ■ ■ 

D i— *■ 

E I— ± 

F h^t 
G ' — ©■ 
H h- • 



300 



-0 — I 



300 



— i — i — i — 



280 



260 



240 



220 



J , , , L 



420 440 460 480 500 
M a (MeV) 



FIG. 8. In this figure we show the values for the mass and width of the a resonance from different papers in the literature, as indicated in the 
plots. In the first two panels, from left to right, the (dark gray) inner strip corresponds to the interval covered by the weighted average whereas 
the (light gray) outer strip is for the mean value, both given in Table IV. In the last panel, the rectangles correspond to the aforementioned 
intervals in the plane. The references are: A [32], B [40], C [120], D [119], E [59], F [118], G [28], and F refers to the NLO 

determination of this work. 



compatible within errors with their weighted average. 6 These 
values are calculated and given in Table IV. 

The ideal situation is that for the threshold parameters ao 
and bo, as can be seen by simple inspection of Fig. 7, or di- 
rectly from the values in Table IV. All values agree within 
errors with their weighted average: 

a° = 0.220 + 0.003 , 
b\Ml = 0.279 + 0.003 . (26) 

The latest NA48/2 Collaboration result [90] is a\ = 0.2220 + 
0.0128 stat + 0.0050 syst + 0.0037 th , in good agreement with 



For a given set of N independent measures .v, with their errors <r,, the 
(weighted) average is given by x = fej^i Xi/&f ) / (Z^Li l/°7j an ^ tne 
standard deviation <x by 1/cr 2 = 2^1] l/trf . 



Eq. (26). For completeness we also report our result at NLO 
for the 1-2 isoscalar scattering length: 

a\ = -0.0424 + 0.0012. (27) 

The last value from K e \ decays of the NA48/2 Collaboration 
[90] is a\ = -0.0432 + 0.0086 stat + 0.0034 syst + 0.0028 th , 
whereas the precise determination of Ref. [35] gives = 
-0.0444 + 0.0010. At this point, it is worth stressing that our 
unitarized amplitudes with the kernels calculated at NLO al- 
low a good reproduction of the low energy behavior (K e 4 data 
and scattering lengths) while keeping the agreement with the 
higher energy data. 

The case of the cr mass and width is not so mild. In Fig. 8 
one can see that the agreement within errors of the different 
values with the weighted average starts at the level of (2 - 
3)cr. At this stage it is then preferable to take the mean of 
the different measures instead of the weighted average. In this 
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way we have: 



= 458 ± 14 MeV , 
r ff /2 = 261 ± 17 MeV . 



(28) 



The resulting error is around 3 times bigger than that for the 
weighted average over the different values considered. The 
different determinations agree within errors with the above re- 
sult, Eq. (28). It can be concluded that our present knowledge 
on the pole position of the cr meson is quite precise, with the 
uncertainty of the order of few tens of MeV, lying in a range 
much narrower than the values nowadays reported in the PDG. 



D. Dependence with M n of the cr meson mass and width 

We can now study the evolution of the cr meson properties 
when the physical pion mass M„ varies (e.g. by changing the 
current quark masses in QCD). This is an interesting problem 
by itself. It is also related to the form factor of the cr me- 
son, Fa-(s), since &s lT l<\M 2 and F lT (0) are proportional by the 
Feynman-Hellman theorem, as discussed below. At LO, the 
only changes produced by varying M 2 are those occurring in- 
side the kernel and the loop function G(s). At NLO, F„ 
varies with M 2 because of Eq. (12), and also the LECs because 
it follows from Eq. (6) that: 



Z i (M^) = / i (M^ phys )-log 



Mi 



M 2 , 

n, pnys 



(29) 



We can consider the subtraction constant a in the function 
G(s) as independent of M n in view of Eq. (22). With the 
above considerations one searches the cr pole position in the 
^-complex plane, s a , for different values of M n , just as in the 
physical pion mass case. The coupling g 2 ^ is also obtained by 
means of the Cauchy theorem. 

Before discussing this evolution, it is useful to make some 
analytical derivations. Let us consider the unitarized nn ampli- 
tude, Eq. (15), as a function of both the Mandelstam variable 
s and the pion mass squared, T(s, M 2 ). In the second Riemann 
sheet it reads: 



t(s,m;) = 



V(s,Ml) 



l+V(s,M 2 )G n (s,Ml) ' 



(30) 



This function has a Laurent series around v expressed in 
Eq. (25). Taking the derivative of T(s, M 2 ) with respect to M 2 
in both sides of Eq. (25), and attending to the double-pole 
terms, one obtains: 

V(M^) = - v ^^ )2 (nv,M„ 2 ) - V( v ,m2) 2 G /7 (v,m2)) 

(31) 

where the dot denotes derivative with respect to M 2 . In the 
previous equation we have taken into account that Eq. (30) 
requires that G//(v) = - ^/V(s,j-) at the pole position s a . 

Analogously, since ^(M 2 .) 2 is minus the residue of the pole 
of the amplitude in the s variable, one gets: 



sliMl) = 



2^2 



V(s a ,M 2 ) 



V'(V,M2) - V(s iT ,M 2 n ) 2 G' n (s iT ,Ml) ' 



(32) 
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FIG. 9. From top to bottom. First (second) panel: Mass (half width) 
of the cr as a function of M n . In the last panel we show the half- width 
as a function of the mass of the a while varying M n . In the figures 
the (red) thinner and (blue) thicker solid lines correspond to the LO 
and NLO results, respectively. In the upper panel the (black) thin 
dot-dashed line represents the two-pion threshold, 2M„. The larger 
circles in the last panel highlight the chiral limit and physical case 
results, whereas the smaller circles represent 25 MeV steps in M„, 
starting at M n = 50 MeV. The dashed, gray lines are the results of 
Ref. [122]. The squares in the first panel correspond to the lattice 
QCD results of Ref. [123], while the rest of points are taken from 
Ref. [124]. 



where the prime denotes a derivative with respect to the s 
variable. One should replace Gii(so-,M 2 ) by G(sa-,M 2 ) (the 
function in the physical Riemann sheet) in Eqs. (31) and (32) 
for the case when the cr pole becomes a bound state. From 
Eqs. (31) and (32), given the knowledge of v and g 2 ^ in the 
physical case, the evolution of the pion pole and the coupling 
with M 2 could be studied directly. We have checked that the 
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M % (MeV) M K (MeV) M K (MeV) 

FIG. 10. From left to right. In the first two panel we show s a in units of the pion mass squared as a function of M n . The second panel shows in 
more detail the region for M n > 300 MeV. In the last panel gj. is depicted in the same units. In the figure the (red) thinner solid lines are for 
the LO results, and the (blue) thicker solid ones correspond to NLO. The solid lines correspond to the real part of the quantity shown, whereas 
the dashed ones represent its imaginary part. We indicate with arrows the points at which s a = AM\ at LO and NLO. 



numerical results are the same as those obtained by looking 
for the pole in the complex plane for different pion masses, as 
explained above. 

The main features of the evolution of the <x meson with M„ 
can be grasped by the inspection of Figs. 9 and 10. In Fig. 9 
we show yfs^ as a function of M„, so that, M a is shown is 
the upper plane, Fo-/2 in the middle one and the plane M ff - 
Tcr/2 in the panel on the bottom. The (red) thinner solid lines 
originate from the LO calculation, V = Vo, and the (blue) 
thicker solid ones from the NLO results, V = V2 + V4, Eq. (18). 
For the physical situation (M„ 140 MeV), we have the case 
just described, that is, the cr meson is seen as a pole in the 
unphysical Riemann sheet. As we increase M n , the imaginary 
part of yfs^ decreases, becoming zero at M n ^ 310 MeV for 
LO and at M„ * 330 MeV for NLO. 7 

In Fig. 10 we show v in units of the pion mass squared in 
the first and second panels from left to right. In the latter the 
scale of the ordinate axis changes and is restricted to values 
slightly slower than AM 2 , so that one can appreciate the evo- 
lution of the real part of s a and distinguish it from the line 
So- - A M 2 (which is difficult to realize from the first panel 
for M„ > 300 MeV). In the last panel we show g 2 ^ in the same 
units for varying M„. For all the panels the solid (dashed) lines 
are for the real (imaginary) part, and the thicker (thinner) lines 
correspond to NLO (LO) results. Notice that both for LO and 
NLO, g 2 T diverges at the point where s a becomes purely real. 
Approaching this point from lower values of M„, Im g^. di- 
verges, whereas, approaching it from higher values of M„ then 
Re gl is the one that diverges. This can be understood from 



7 At this point another pole (not shown in the figures) starts to appear be- 
low the cr one. This is due to the appearance of two real solutions for the 
equation 1 + V(s)G(s) = 0, since the imaginary part of s& is zero in this 
region. There is no need to consider further this pole since, irrespectively 
of whether it lies in the same Riemann sheet than the higher pole, the ef- 
fects of the latter overwhelmingly dominate over those of the former. For 
smaller M n , since the solutions are not real, the tr corresponds to two com- 
plex conjugated values. 



the behavior of the derivative of v, that is not defined pre- 
cisely at this point, and in view of Eq. (31), where it is seen 
that So- oc g 2 a . 

For even larger values of M„ (M„ 370 MeV at LO and 
M n a 480 MeV at NLO), s lT osculates the 2n threshold, while 
standing below it, and changes from the unphysical Riemann 
sheet to the physical one, becoming a bound state. Since s a 
AM 2 close to this point, the binding energy is small, and so is 
the coupling, becoming exactly zero when s lT - AM 2 . These 
points are indicated with arrows in Fig. 10. This behavior can 
be shown analytically. From Eq. (32), one deduces that for 
Sa - AMI 

gl = -77 6AnM n ^|v-4M2| , (33) 

with j] — +1 for the unphysical Riemann sheet (at the left 
of this point) and r\ = — 1 for the physical Riemann sheet 
(at the right). Therefore, g 2 a — for v = AM 2 , as indi- 
cated by the arrows in the rightmost panel of Fig. 10. How- 
ever, it is worth noticing that from Eq. (33) it follows that 
gl/ yj\ v/4 - Ml\ = g^/lpo-l is finite. On the other hand, the 
fact that the pole changes from one Riemann sheet to the other 
in a continuous way can be understood in terms of Eqs. (23) 
and (24). The difference between the G{s) function calculated 
in the two Riemann sheets is given by a piece proportional to 
o-(so-) = yl -AM 2 /s lT that vanishes for s a = AM 2 . At this 
point, where the cr is a zero bound state, one also has an infi- 
nite value for the scattering length. 

The mere existence of this critical point can be examined 
analytically. For s = AM 2 , the function G{s) can be written as: 



GO 



AM 2 ) 



a + log 
16^2 



log 



16tt 2 



(34) 



with n 2 = e~ a fi 2 a new scale. If we concentrate on the simpler 
case of LO, V(4M 2 ) = 1M 2 I2F 2 , the equation for finding a 
pole at s = AM 2 ., V' 1 + G = 0, can be cast as f(x) = 0, with 



f(x) — 1 + ax log x , 



(35) 
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FIG. 11. Representation of the function f(x), Eq. (35), for two values 
of a, a > a* (upper line) and a < a* (bottom line). The variable x is 
defined in Eq. (36). 



where 



and 



lnll{32n 2 F 2 n ) > 



(36) 



(37) 



Since a > a zero of the f(x) function is only possible for 
< x < 1. Actually two zeros of this function exists if the 
value of the function at its minimum xq — e~ l is negative (see 
Fig. 1 1). This condition in terms of the variable a requires that 
the latter is smaller than the critical value a*, 



-1 +log 



7H 2 
32n 2 Fl 



(38) 



If this is the case there is a zero for < x < Xq and another 
one for xo < x < 1 . For our value of the renormalization scale, 
ju = 770 MeV, a* =* -0.6, so that the fitted value a = -1.4 
given in Table I is much smaller than a*. We also have that 
our value for xq corresponds to M„ 900 MeV, then a pole 
with So- = AM 1 - exists for < M n < 900 MeV. The solution of 
Eq. (35) for the value of a fitted gives that this pole is located 
at M„ a; 370 MeV, as stated above and indicated by the left 
most arrow in the panels of Fig. 10. 

For the NLO case, the situation becomes somewhat more 
involved, and the function f(x) is now: 

f{x) = 1 + a(x)xlogx\l + a(x)J3(x)x) , (39) 

where a(x) is defined as in Eq. (37), but at NLO one has to 
take into account its implicit dependence onioc M 2 through 
F„. On the other hand, /3(x) is defined as 

40 2 

B(x) — L log x , 

P ' 147 7 s 

L = S + 2?-!5 + g5 + ! + ^lpg*„ (40) 

where = li(M 2 phys ) corresponds to the LECs calculated at 
the physical pion mass and 



For the values collected in Table I we find that s a = 4M% 
for M„ =i 480 MeV. Nevertheless, this value is quite sensi- 
tive to the LECs, and it should be taken merely as indicative 
(for some values of the LECs not far from the fitted ones the 
change from virtual to bound state does not occur at all). This 
sensitivity is illustrated by the error band in Figs. 9 and 10. 

In Fig. 9 our results on the pion mass dependence of the 
cr pole position, partially presented in Ref. [121], are com- 
pared with other works. The (gray) dashed line, denoted by 
IAM, gives the results of Ref. [122] in the framework of the 
IAM. The points shown come from the lattice QCD stud- 
ies of Refs. [123, 124]. Interestingly, we find a remarkably 
good agreement with the curve from the IAM results [122] 
for M„ < 400 MeV. As stated by the authors, the point where 
s a = 4M^, and thus the cr meson becomes a bound state, is 
M n a 460 MeV when they employed the NLO ChPT ampli- 
tudes [125], whereas M„ 290-350 MeV when the N 2 LO 
ChPT amplitudes were used. We show in Fig. 9 the curves of 
Ref. [122] corresponding to this latter case. 

A lattice QCD search of light scalar tetraquarks with J PC = 
++ (we focus here on the 7 = results) is performed in 
Ref. [124]. Along with the lowest 7r(p)7r(-p) scattering state, 
an additional lighter state is found. For the dynamical sim- 
ulations of Ref. [124] the former state is denoted in Fig. 9 
with « = 1 (green filled circles) and the latter one with n — 2 
(pink filled triangles). For the quenched simulations we use 
the (green) empty circles and the (purple) empty triangles, in 
the same order as before. 8 The points with n = 1 and 2 over- 
lap at each pion mass, and the quantitative agreement with our 
curves is satisfactory. However, both our curves and the lattice 
QCD results of Ref. [124] do not agree with most of the points 
of the lattice QCD calculation of Ref. [123] and, in addition, 
the tendency of the points is qualitatively different to that for 
our results and those of Ref. [124]. 

For larger values of M„ we obtain values for the cr meson 
mass, both at LO and NLO, that remain below but always 
close to the 7T7t threshold, in agreement with the lattice QCD 
results of Ref. [124]. Note that this is not the case for the IAM 
calculation of Ref. [122] for M„ > 400 MeV. The fact that 
the cr meson follows so closely the threshold for higher val- 
ues of M„, both according to our calculation and to the lattice 
QCD calculation of Ref. [124], clearly indicates that for such 
masses it is dynamically generated from the 7T7t interactions. 
We elaborate further on the nature of the cr resonance below. 
However, one should keep in mind that the cr meson becomes 
an anti-bound or virtual state between those pion masses in 
which it has zero width and has not crossed to the physical 
Riemann sheet yet. In the bound state case, an additional state 
appears in the energy levels spectrum in the box, whereas an 
anti-bound state does not. In order to discern the latter situa- 
tion one should look at other computable quantities, such as 
the sign of the / = S -wave nn scattering length. 



X P = M lphyJl4 ■ 



(41) 



' However, we must also point out that the lattice QCD simulations are per- 
formed for each pion mass at a single volume and lattice spacing, so the 
continuum and infinite volume values of the o" meson mass in the bound 
state case may differ from those values. 
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FIG. 12. Kinematics of the nn scattering process in the presence of 
a scalar source, n{pi)n(p2)H{q) — > n{pi)n(p^). Pions correspond to 
the solid lines and the scalar source to the wavy one. The gray blob 
indicates the interactions involved. 



It is also interesting to study the chiral limit, M n — > 0. As 
can be seen in Fig. 10, Sa-jM 2 — » oo, because remains finite 
in this limit. Indeed, the values calculated for v near the chiral 
limit behave as (for M„ < 150 MeV), 



Sa-iMl) = Sa-x +aMl+bM 2 n log 



Mi 



M 2 . 

7r,phys 



(42) 



with the values of the cr pole position in the chiral limit 



given by 



453 - /282 MeV (LO) and 



402 - /263 MeV (NLO), see Fig. 9. 



IV. THE SCALAR FORM FACTOR OF THE cr MESON 

We turn now our attention to the calculation of the scalar 
form factor of the cr meson, that is, the interaction of the a res- 
onance with a scalar source (denoted in the following by H). 
As an intermediate step we calculate first the scattering of two 
pions in the presence of a scalar source, from which we extract 
the scalar form factor of the cr. This can be done because the 
cr originates as a pole in the interaction of a scalar isoscalar 
pair of pions, as discussed in Sec. III. We start by considering 
in Sec. IV A the kinematics of the ttttH — > nn reaction, which 
is somewhat more complicated than the standard kinematics 
of a two-body reaction. In Sec. IV B, we discuss the one-loop 
calculation of the amplitude nnH — > nn from the chiral La- 
grangians of Sec. II. In terms of this amplitude one can de- 
rive the scalar form factor of the cr meson, as performed in 
Sec. IV C. This is accomplished by taking into account pion 
rescattering, similarly as done above for nn scattering, with 
some modifications that are carefully examined. 



A. Kinematics 

We are interested in pion-pion scattering with a scalar 
source, n(p\)+n(p2)+H(q) — > n(pi,)+n(p4), Fig. 12. The over- 
all center-of-mass frame, CM, is the same as the rest frame of 
the final pions, while that corresponding to the initial ones is 




FIG. 13. The unit three-momenta in terms of the polar and azimuthal 
angles. 



denoted by CMB. Due to the presence of the scalar source 
CMB does not coincide with CM. In the CM one has 



(43) 



and 



Pi + P2 = -q 

q = (q°,q). 



(44) 



We denote by s and s' the invariant masses squared for the 
final and initial pions, in order. At the end of the calculation, 
the limit s, s' — > .v is taken. It follows that 



(Pi + Pif = s' = (^-q°) -q 2 = S + q 2 -2q°^, (45) 
and then, 

s — s' + q 2 



'I 



2yfs ' 



(s + s'-</ 2 ) 2 

q = -, s 

4s 



(46) 
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Analogously, one has in CMB: 



PI= 


I z / 


P2 = 


{ L 1 




s' 

--Ml, 
4 


«° = 


s — s' — q 2 




q* = 


(s + s' - q 2 ) 


4s' 



(47) 



In the following quantities with a tilde are expressed in CMB. 
Notice that p is the three-momentum of the first pion in CMB, 
while p refers to the three-momentum of the third pion in CM. 

The final (initial) two-pion states are projected into S -wave 
in CM (CMB) because the cr resonance is defined as a pole 
in the second Riemann sheet of the mi isoscalar S -wave in 
CM (CMB). The unit three-momenta (indicated with a hat) are 
given in terms of the polar and azimuthal angles (see Fig. 13) 
as: 

p = (sin 8 cos <p, sin 6 sin <f>, cos 6) , 
p = (sin § cos (f>, sin 9 sin 0, cos , 
q = (0,0,1), (48) 

where we have chosen the z-axis to be the direction pointed by 
q. We now work out the Lorentz transformation from CMB to 
CM: 



(Pi +P2)cMB = (V?,5) , 

(pi + pi)cm = (Vs-?°.-q) = ( 



S + S - q- 
2^ 



-q).(49) 



The transformation reads: 



s + s — q 



2^~s 



yV7, 

q = -y Vs'v 



It follows then that y =11 Vl - v 2 and v are 



y 



s + s - q 
q 2 yfs 



q 



y V7 S + s' - q 2 
We further define the four-momenta S and A given by 

£ = (Pi +P2) , 

A = (pi- p 2 ) . 



(50) 

(51) 
(52) 



In the CM 



I = (VW,~q). 



(53) 



(54) 



The momentum transfer A has a simple expression in CMB 
where it is given by A = (0, 2p). We then perform its Lorentz 
transformation to CM, with the result 



A u = -2 



q p 



2p + 2(y 



l)[pl(p-q)q. 



(55) 



The problem has six independent Lorentz invariant kine- 
matical variables. 9 We define, in analogy with two body scat- 
tering, the following six alike Mandelstam variables, 



(P3 


+ P4?, 


Q>1 


+ Plf, 


(Pi 


-Plf, 


(P2 


-PA) 1 , 


(PI 


-Pa) 2 , 


(P2 


-P3) 2 . 



These variables fulfill the relationship 



t + u 



mi 



(56) 



(57) 



+■ t +u = q - 

which is the analogous one to s+t+u = 4M 2 . valid for two-pion 
scattering, Eq. (13). Though q 2 and the variables in Eq. (56) 
are not independent because of Eq. (57), it is convenient to 
write the different amplitudes nnH — > nit in terms of all of 
them, given the symmetries present in the calculation. 

In virtue of the previously worked Lorentz transformation, 
Eq. (52), we have the four-momenta properly defined in CM 
in terms of the key variables s, s', q 2 and the polar and az- 
imuthal angles in the two-pion center of mass frames (the 
Lorentz invariants only depend on the difference between the 
azimuthal angles, see Eq. (59) below). It is convenient to ex- 
press p\ = (S + A)/2 and p 2 = (£ - A)/2, with S and A given 
in CM by Eqs. (54) and (55). In terms of this set of variables, 
the Lorentz invariants of Eq. (56) are given by 

t = 2M 2 n -2(a + A + B + C) , 
f = 2M; - 2{a-A-B + C) , 

u = 2M 2 „ -2(a+A-B-C) , 

u = 2M 2 -2(a-A+B-Q 



(58) 



where 



A = -A 
2 



pIa = ~lq 



r- ~ 
[p|— — COS0 , 

V7 



B = --t 



P = +2WHPI cos( 



c 



1 



A - p = -|p||p| 



p p + — 



2yfsV? 



(59) 



■ COS0COS0 



9 One the five four-momenta involved in the reaction is fixed by energy- 
momentum conservation. From the other four ones we can construct 6 in- 
dependent scalar products. Notice that p\= p\ = p\ = M% and that g 2 can 
be derived from Eq. (57) in terms of other Lorentz invariants. 



In the previous equation the five kinematical variables, s, 
s', q 2 , cos 9, cos 9 are used together with the scalar product 

p • p = sin 9 sin 0cos(</> — 0) + cos 9 cos . (60) 

In terms of the variables in Eq. (56) one can express the in- 
verses of several pion propagators that appear in many Feyn- 
man diagrams that contain the scalar source attached to an 
external pion leg, cf. diagram (a. 2) of Fig. 14. It results: 



D, 


= (q + pi) 2 


-M\ 


— s 


+ t' 


+ u' 


-AMI 


D 2 


= (q + Pif 


-Ml 


— s 


+ t 


+ It 


-AMI 


D 3 


= (q- Pi) 1 


-Ml 


= s' 




+ u 


-AMI 


D 4 


= (q- p\f 


-Ml 


r 

= S 


+ t 


+ u 


-AMI 



Because of Eqs. (59), the angular dependence of these inverse 
propagators is rather simple: D\ and D2 depend just on cos 9, 
while D3 and D4 do on cos 9. The propagating pion can be- 
come on-shell for certain angles, giving rise to a pole in the 
propagators. These poles, when the S -wave angular projec- 
tions are performed, result in logarithmic divergences. In par- 
ticular, there is always a pole for q 2 — > 0. We treat this issue 
later on. 



B. The nn H — > mt scattering amplitude 

To determine the Feynman diagrams required for the nn 
scattering in the presence of a scalar source up to 0{p A ) in 
ChPT it is useful to have in mind first those diagrams of 
plain nn scattering in Sec. Ill, Fig. 2. Now, one external scalar 
source has to be added in all the possible ways to those dia- 
grams. As deduced from the Lagrangians X.2 and £4, Eq. (2), 
the scalar source can couple to any even number of pions. 
In Fig. 14 we show the diagrams that must be calculated at 
the one-loop level, where the external scalar source is indi- 
cated by a wiggly line. The LO diagrams correspond to (a.l) 
and (a. 2). 10 Diagrams (a. 2), (e.l) and (f.l) can be handled to- 
gether because their sum correspond to taking the full pion 
propagator in between the external source and the four-pion 
vertex, Eqs. (10) and (11). In addition, all the diagrams on the 
bottom line of Fig. 14, namely, (e.l)-(O), correspond to the 
wave function renormalization of the LO ones. Both issues are 
derived to NLO from the pion self-energy diagrams, Fig. 1, 
Eq. (1 1). Once the the renormalization of the pion propagator 
and the the wave function renormalization are taken into ac- 
count, as well as the rest of diagrams diagrams in Fig. 14, one 
has the basic topologies shown in Fig. 15. 

Compared with nn scattering the presence of the c-number 
external scalar source H complicates considerably the simple 
expressions for the former [7]. The calculation for each of the 
diagrams in Fig. 15 is given in Appendix B. Specifically, we 



Of course, the scalar source can be attached to any of the pion legs but for 
conciseness we draw explicitly the attachment to only one. This should be 
understood in the following. 
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(c.l) (c.2) (c.3) 




(d.l) (d.2) (d.3) 




(e.l) (6.2) (6.3) (f.l) (f.2) (f.3) 



FIG. 14. Feynman diagrams for the nn scattering amplitude in the 
presence of a scalar source, nn H — > nn, at one-loop order in ChPT. 




1 II III IV V 




VI VII VIII IX X 




XI XII XIII 



FIG. 15. Final set of Feynman diagrams for the nn scattering in the 
presence of a scalar source, nn H — > nn, at 0(p 4 ) in ChPT omitting 
the pion propagator dressing and wave function renormalization of 
the leading order diagrams in Fig. 14. 

calculate the processes n Q (p 1 )n°(p2)H(q) — > n°(p3)n°(p4) and 
n°(p 1 )n Q (p2)H(q) — > n + (pi)n~(pn,), with the former denoted 
by T„„ and the latter by T„ c . These two processes are con- 
sidered in order to isolate the pion pairs with definite isospin 
(/) by taking the appropriate linear combinations. The stan- 
dard decomposition of the n°n° and n + n~ states in two-pion 
isospin definite states, \nn{II{)}, being I3 the third-component 
of isospin, is 

kV>= ^||^(20)>- ^|7T7r(00)), 

\n + n-} = - ^||;nr(20)> - ^|7T7r(10)> - ^\nn(00)) , 



16 



W(s, q 2 ) 





■V(s')G(s')S l + V(s)G(s) 



FIG. 16. External scalar source coupled to a double a pole in the 
mt H — > mi process. The <r pole is originated by the resummation of 
pion re-scattering, as indicated in the left diagram by the iteration of 
the unitarity two-point function. 



where we have taken into account that \n + ) = -\tt; 1=11^ = 
-1), as follows from the definition of the n + field, Eq. (5). 
Because of isospin conservation (the scalar source H(q) is 
isoscalar), the Wigner-Eckart theorem implies 

{n a n°\S\n n° s) = + |<7r7r(20)|^|7T7r(20)//) 
+ ^{7m(QQ)\S\nn(QQ)H) , 

{n + n-\S\n°n° s) = - -(7m(20)\S\nn(20)H) 

+ -(im(00)\S\nn(00)H) , (62) 

with S the S -matrix. From this equation we can isolate the 
purely 1 = matrix element, &{s, s' , q 2 , 6, 6, <f>, 0), corre- 
sponding to 

M(s, s', q 2 , 0, 8, <f>, 4>) = {7T7T(00)\S\7m(00)H) . (63) 

From Eq. (62), we have: 

Jl(s, s', q 2 , 8, 8, cf>, 0) = {n Q n°\S\n K° s) + 2<7r + 7r|.S|7rV s) 
= T m + 2T m . . (64) 

We are interested in this matrix element because the cr is 
isoscalar. 

The cr is an S -wave resonance so that it is also required 
the S -wave angular projection of the initial and final isoscalar 
pion pairs. This is straightforward for the final pions because 
the CM coincides with its own rest frame, with the result: 

\tt7t- 00)34 = ^ jdp \n(p 3 M P4 )(00)) . (65) 

Regarding the initial pair of pions, its state is defined in CMB 
analogously as in the previous expression. One has still to per- 
form the Lorentz boost to the CM frame so that 

Itttt; 00> 12 = f/(v)^ j 4\7T(^f,p)7T(^f,-p)(00)), 

(66) 

where U (v) is the Lorentz boost operator from CMB to CM, 
with the velocity v given in Eq. (52). When acting on the pion 
states (which have zero spin) the only effect is the transfor- 
mation of the four-momenta from CMB to CM. Then, we can 
also write Eq. (66) as 

\7t7t- 00} n = ^ J 4^(/>i(p)MP2(p))(00)> , (67) 



where p\ and p2 are written in terms of the four- momenta in 
CMB. From Eq. (53), pi = (I + A)/2, p 2 = (2 - A)/2 with 
Z and A given in Eqs. (54) and (55), in order, as a function of 
the CMB kinematical variables. 

Employing the states projected in 5 -wave, Eqs. (65) and 
(67), we are then ready to calculate the required matrix ele- 
ment in ChPT, if{s, s' , q 2 ): 



<p(s,s',q ) 



" 32tt2 / 



d 2 Q 3{(s, s',q 2 ,6, 8, 



(68) 



Note that the extra factor 1 /2 in Eq. (68) arises because of the 
unitary normalization, as explained after Eq. (14). In the last 
equation, the double solid angle integration is 



/m r*Ln r*Ln 

d 2 Cl= ( dcos6» I dcosg ( d0 I d0 . (69) 
J-i J-i Jo Jo 

One linear combination of azimuthal angles, <p an d ^> is a 
spare variable, and then one integration in Eq. (69) is triv- 
ial. This is so because they appear just through the expression 
cos(0 - 0), as explained above, Eq. (59). In fact, for any peri- 
odic angular function, f{y) = f{y + 2n), one has: 

J fin fin fin 
dcf, d0/(0-0) = 27T d r /( T ). (70) 
o Jo Jo 



C. Scalar form factor 

Once the perturbative amplitude for the process nn H — > nn 
is calculated, we proceed by taking into account pion rescat- 
tering, similarly as was done for mt — > nn, see Eqs. (15)- 
(18). The resulting amplitude is denoted by T$(s, s',q 2 ), and 
following the same unitarization method as in Sec. Ill from 
Refs. [19, 24], it can be written as: 



T s (s, s',q 2 ) = 



W(s, s',q 2 ) 



(1 + V(s)G(s))(l + V(s')G(s')) 



(71) 



This is the analog to Eq. (15) but now for the process 
jxnH — > nn, with the new kernel W(s, s',q 2 ) instead of V{s) 
in Eq. (15). It is important to stress the presence of two fac- 
tors 1 + VG in the denominator of Eq. (71). This is so because 
in 7T7r H — > 7T7T the presence of the scalar source H{q) makes 
necessary to resum the unitarity loops corresponding to both 
final and initial state interactions. 

The kernel W(s, s', q 2 ) is obtained in a chiral expansion by 
matching Eq. (71) order by order with its perturbative calcu- 
lation. The chiral expansion of the kernel is 



W = W 2 + W 4 + 0{p b ) , 



(72) 



where we omit the dependence on the arguments s, s' and q 2 
for easy reading. The subscripts in Eq. (72) refer to the chiral 
order. Then, the amplitude Eq.(71) is expanded, as it was done 
in Eq. (17), so that one has: 

T s (s, s',q 2 )=W 2 + W 4 

-W 2 V 2 (s)G(s) - W 2 V 2 (s')G(s') + 0(p 6 ) 
=<p 2 + <P4 + 0(p 6 ), (73) 
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where <p„(s, s',q 2 ) is the 0(p") contribution to ip(s, s', q 2 ) de- 
fined in Eq. (68). The kernels W„(s, s', q 2 ) are determined by 
matching the above expressions order by order, so that: 

W 2 (s, s',q 2 ) = <p 2 (s, s',q 2 ) 

W 4 (s, s', q 2 ) = <p 4 (s, s', q 2 ) + <p 2 &(j)G(j) + V 2 &(S')G(S') , 

(74) 

where it was used that V2(s) = %2(s), Eq. (18). 

The form factor of the <x meson, FAq 2 ), can now be ex- 
tracted from Ts(s, s', q 2 ), employing W — W2+W4 in Eq. (71). 
For that one has to isolate the double cr pole present in 
Ts(s, s',q 2 ), as drawn on the right-hand side of Fig. 16. The 
double cr-pole contribution can be written as: 11 

FAq 2 ) gl ]im Wjs, s',q 2 ) 

(s - softs' - s^) >S->* (1 + V(s)G(s)) (1 + V(s')G(s')) 

(75) 



Expanding the r.h.s. of the above equation around s, s' 
So-, and equating the double-pole term, the result is: 



FAq 2 ) 



Visa-} 



■W(v,V,<7 )• 



(76) 



In determining the kernels W n (s,s',q 2 ), we have followed 
the master guidelines of pure jttt scattering procedure to take 
into account the rescattering of the pions given in Sec. Ill A. 
However, some modifications are needed in our case because 
of the presence of the pion propagators in the external pion 
legs attached to a scalar source, see Fig. 15. Let us focus, for 
clearness, in the LO amplitudes (pi{s, s' , q 2 ), corresponding to 
the diagrams I and II in Fig. 15 (the amplitudes are given in 
Appendix B). Before the angular projection in Eq. (68), one 
has 



M 2 (s, s',q 2 ,0, S,( 



2B 

'n 



2 X 



Ml/2 



Di 



where the subscript in J[ refers to the chiral order, = s' 
and 53 4 = s, and the D, are the inverse of the pion propa- 
gators given in Eq. (61). These contributions proportional to 
the propagators stem from the piece of diagram II in which 
the on-shell part of the An vertex is retained, so that the pion 
propagator is not cancelled out by an off-shell part from the An 
vertex (cf. Ref. [13]). Considering, for conciseness, the case 
s — s' (the one interesting for the cr scalar form factor for 
which s = s' = So-), these propagators can be written as: 



1 



1 



f -2|p||q| cosfj 



(78) 



and similarly for the other D,. It should be noted that for cer- 
tain values of q 2 and s, these propagators can have a pole in 
the variable cos 6. In particular, for q 2 — > this is always the 




FIG. 17. String of unitary diagrams to be resummed when the scalar 
source is attached to an external leg. On the right, the resulting single 
cr pole contribution is depicted. 



case. Upon angular integration, this contribution gives rise to 
an imaginary part that diverges as 1/ y/\q 2 \ for q 2 — > 0~. As 
shown below, this limit is the one that matters in order to cal- 
culate the quadratic scalar radius of the cr, but this divergence 
would lead to an undetermined value for it. This fact is not ac- 
ceptable and indicates a deficiency in the procedure followed 
up to now. 

Let us clarify this important technical point and the 
way it can be solved. The term of the amplitude 
^2(s,s',q 2 ,6,6,(p,^)) in Eq. (77) that is proportional to the 
pion propagators, that we denote by ^2,prop, 12 can be written 

as: 



prop 



^ AB 



(79) 



where we have taken into account that £2 (*) = {s-M 2 j2)jF 2 n . 
Once J?l2,piop is projected in the S -waves for the initial and fi- 
nal pion pairs, Eq. (68), we end with the contribution W2, prop 
to the kernel W{s, s, q 2 ) in Eq. (71). Keeping in this resum- 
mation only terms up to one-loop, and hence proportional to 
G(s), one obtains 7 , S 2, P iop given by 



7s2,prop = -2£>(s)G(s)VK2,prop(s, <? ) , 



(80) 



(77) where the expansion 



1 



2 = Z (_1) " ( " +1)( ^ )G(s)) "' (81) 



n=0 



is employed. 13 

However, the result of the one loop calculation in ChPT 
of the diagram VIII, once properly projected in isospin and 
S -waves as discussed above, gives half of the amplitude in 
Eq. (80). Whence Eq. (71) is double counting this kind of 
terms at the one-loop level. Analogously, it can be seen in the 
n-loop terms of the resummation that the contribution of the 
kernel proportional to W2, pm p is counted n + 1 times, Eq. (81). 
This is so because we are missing the proper combinatoric 
factors as an on-shell factorization scheme for unitarizing 
is employed. Thus, instead of resumming these terms with 
1/(1+ VG) 2 , they should be resummed with just 1/(1+ VG) in 



12 Recall that we are interested in the s = s' case. 
1 1 Because of invariance under temporal inversion the amplitudes for nn — * cr 13 Recall that V(s) in Eq. (71) is ^2(*) because we are unitarizing a one-loop 
and cr — > nn are equal. ChPT calculation for nn H — > nn. 
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order to give the proper diagram counting. Notice that in this 
case they do not contribute to the double-pole term needed for 
Fo-(q 2 ), as can be seen from Eqs. (75) and Eq. (76). This is 
also shown schematically in Fig. 17. Had we considered an 
integral equation for the resummation procedure instead, this 
kind of contributions would be integrated giving terms pro- 
portional to the three-point function Co{s, s', q 2 ), in which the 
scalar source interacts with intermediate pions, like the terms 
appearing in the diagram X of Fig. 15. This is not a short- 
coming of our approach, because this kind of diagrams are 
properly included when the kernel W(s, s', q 2 ) is calculated at 
higher orders in the chiral counting, as can be seen in Fig. 15. 
E.g. at the one-loop level calculation of W(s, s',q 2 ) one has 
the diagram X of Fig. 15, that arises from iterating once the 
pion-propagator contributions at tree level. 

From the previous discussion we remove the terms of the 
amplitudes with the external scalar source coupled to initial 
or final asymptotic pion legs from the kernel W(s, s',q 2 ) in 
Eq. (71), as they do not contribute to the scalar form fac- 
tor of the cr. The latter requires the coupling of the exter- 
nal scalar source to intermediate pions and vertices. Now the 
question arises of how to remove properly the terms arising 
from the Feynman diagrams with the scalar source attached 
to a pion propagator in an external pion leg. We cannot sim- 
ply drop these diagrams because the pion propagator between 
the source coupling and a pure pionic vertex in an external 
pion leg may be cancelled by off-shell terms from the mr in- 
teraction vertex [13]. Indeed, such contributions are required 
in order to have results independent of pion field redefinitions 
that mix diagrams with different number of pion propagators. 
Rather, a procedure based on the full on-shell amplitude cal- 
culated in ChPT up to some order, which is independent of the 
former redefinitions, must be given. 

Let us consider the general case, and write these contribu- 
tions as: 



f(x,y) 
x — xq 



(82) 



where x = cosf? and xo = g ,2 /(4|p||q|). 14 Here we have col- 
lected in y the rest of the variables. In order to subtract the 
pure pole contribution in Eq. (82) we subtract from the nu- 
merator above the residue of the pole, 



f(x,y)- f(x ,y) 



x- x 



(83) 



In the LO case, in view of Eq. (77), this amounts to removing 
the whole term proportional to the propagator, since it just de- 
pends on s (or s') and not on 9 (or 9), that is, df(x,y)/dx = 0. 
This subtraction procedure is independent of pion field redef- 
inition because in f(xo,y) all the pion lines are put on-shell so 
it cannot contain any off-shell remainder that could be coun- 
terbalanced by other off-shell parts coming from other ver- 
tices, and giving rise to possible pion field redefinition depen- 
dences. 



14 We are considering again the case in which the scalar source is attached to 
n(pi ), since the argument for the the other cases is analogous. 



With this procedure we are then ready to calculate F^iq 2 ). 
For that we define the new amplitude &(s, s', q 2 , 9, 9, 0, 0) ob- 
tained from the original Jl{s, s' , q 2 , 9, 9, 0, 0), Eq. (64), by re- 
moving the contributions with the scalar source attached to 
an external pion leg, following the procedure in Eq. (83). In 
terms of the former we calculate its angular projection as in 
Eq. (68), obtaining the new amplitude <&(s, s', q 2 ): 



®(s, s',q 2 ) 



32n 2 J 



d 2 Cl&(s, s',q 2 ,9, 0,0,0) 



(84) 



Then, the final expression for the interaction kernel, that we 
now denote by 'Wis, s',q 2 ), is (cf. Eq. (74)) 



'W ^ f W 2 + f W A , 
f W 2 = $2 , 

W 4 = $4 + <D 2 &(j)G(j) + 4> 2 fr(s')G(s') , 

with the subscripts indicating the chiral order as usual. 
The scalar form factor of the cr is finally given by 



F<r(q 2 ) = 77} ^To WCv, So-, q 2 ) ■ 



(85) 



(86) 



For definiteness let us explicitly give the expressions at LO 
and NLO for F lT (q 2 ) from the previous equation: 



(s L0 ) 2 

^V)=^T^2(V,V^ 2 ) 



V 2 (V) 2 



With 



(v 2 (scr) + v 4 (v)r 

x(<W 2 (s IT ,s a ,q 2 )+'W 4 ( 2 )) ■ (87) 

V 2 C0 



(^ LO f = lim (^ L0 - s) ^ 



V 2 (s) + V 4 (s) 



(V 2 (s) + V 4 (s))G(s) 



(88) 



where s]f and s^ h0 are the cr pole positions at LO and NLO, 
respectively, given in Table III, and for V 2 and V 4 see Eq. (18). 

One technical detail is in order. The cr resonance is a pole 
in the second Riemann sheet of nn scattering for the physical 
pion mass. As we have seen in Sec. HID when increasing the 
pion mass above some value the cr meson becomes a bound 
state and moves into the first Riemann sheet (the correspond- 
ing pion mass value is indicated by the arrows in Fig. 7). Then, 
Eq. (86) has to be understood in the same Riemann sheet as the 
cr pole happens. This requires the evaluation of "W(s, s', q 2 ) in 
Eq. (85) either in the first or second Riemann, according to the 
value taken for the pion mass. 15 

We now discuss the analytical continuation of the loop 
function Cq(s, s, q 2 ) to the second Riemann sheet (we take 



15 This qualification is only relevant for "W^is, s, q 1 ). 
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FIG. 18. The normalized scalar form factor of the cr meson calcu- 
lated at NLO for the physical case with yfs^ given in Table III. 
The range in q 2 extends from q 2 -0.6 GeV 2 s -30M 2 up to 
q 2 a 0.08 GeV 2 * AM 2 . 



from the beginning in the present discussion that s' = s), 
where it is denoted by Co ; //(s, s, q 2 ). The function Co{s, s', q 2 ) 
corresponds to the three-point one-loop function of diagram 
X in Fig. 15 and its calculation is discussed in Appendix A. 
In order to proceed with the analytical continuation we first 
evaluate the difference 

AC(s, q 2 ) = Cq(s + ie, s + ie, q 2 ) - Cq(s - ie, s - ie, q 2 ) 

(89) 

for s and q 2 real and q 2 < AM 2 } 6 The second Riemann sheet 
in 7T7T scattering is reached by crossing the real s-axis above 
threshold, s > AM 2 , and so we have to consider Eq. (89) 
for the same values of s. It turns out that a cut in s extends 
2Ml 



for 



s > 



M n ^AMl 



■q 2 ee 



s K for which AC(s, q ) 



is non-zero (the same expression for the cut also occurs for 
s < 2M 2 - M„ ^AM 2 - q 2 ). When q 2 -> + (this limits gives 
the same value for the quadratic scalar radius as q 2 — > 0~) to 
cross the real axis for s > AM 2 implies to consider AC(s, q 2 ) 
given by the mentioned cut for Co, s > s IC , correspond- 
ing to AfcCrj in Eq. (A17). Once this discontinuity is evalu- 
ated we continue it analytically in s and q 2 and subtract it 
to Cq(s, s, q 2 ) (calculated in the first Riemann sheet), as done 
above to determine Gu(s), Eq. (23). It results, 



Cojiis, q 2 ) = C (s, q 2 ) - AC(s, q 2 ) . 



(90) 



Notice that for calculating Eq. (85), it is not necessary 
to use Gji(s) when the cr pole remains in the second Riemann 
sheet. This is due to the fact that $4 contains the two-point 
one-loop function Bo(s), evaluated in Appendix A, so that the 
discontinuity when crossing the unitarity cut above threshold 
cancels mutually between these two functions. 



16 For q- > AM 2 the opening of the 2n production process introduces addi- 
tional complications that we skip now since we are mostly interested to 
values of q- around zero, used below to calculate the quadratic scalar ra- 
dius of the cr resonance. The whole region q 2 < 4M 2 is of interest and 
considered by us as well. 



We show in Fig. 18 the modulus squared of /v(<? 2 ) nor- 
malized to Fo-(O) for q 2 < AM 2 calculated at NLO with the 
physical value of M„. We observe a monotonous increasing 
function with q 2 . The LO result is just a constant because <t>2 
is so and is not shown in the figure (it would be just 1). 



V. QUADRATIC SCALAR RADIUS OF THE a MESON 
AND THE FEYNMAN-HELLMAN THEOREM 

The quadratic scalar radius of the cr resonance, (r 2 )^, is 
related to the scalar form factor of the cr by a Taylor expansion 
around q 2 = 0, 



FAql = FA0) + 



dFAq 2 ) 



dq 2 

= 7v(0)(i + |v>: 



q 2 =0 



(91) 



where the ellipsis indicate higher powers of q 2 in the Taylor 
expansion. In this way, 



6 OFAq 2 ) 



FA0) dq 2 



(92) 



q 2 =0 



Notice that, since the form factor reduces to a constant (in- 
dependent of q 2 ) at LO, we find that (r 2 )^" = for this case, 
similarly as the case for the quadratic scalar radius of the pion 
[7] within standard ChPT. Whence, the quadratic scalar ra- 
dius must be calculated at least at NLO. Before discussing the 
results for the physical pion mass case, we study the depen- 
dence of (r 2 )^ with the pion mass. We show the square root 
of the quadratic scalar radius of the cr, y/{r 2 )f, in the upper 
panel of Fig. 19 as a function of M„, with its real part given by 
the (blue) solid line and its imaginary part by the (red) dashed 
line. It diverges in the chiral limit (M„ = 0) and where the cr 
pole coincides with the two-pion threshold (indicated by the 
rightmost arrow in Fig. 10). The latter point corresponds to a 
zero energy bound state and as such it must have infinite size, 
as dictated by elementary quantum mechanics. On the other 
hand, in the chiral limit (r 2 )^ also diverges as log M„, simi- 
larly as the quadratic scalar or vector radius of the pion [7], 
because the infinite size of the pion cloud around the bosons. 
It is relevant to note that the imaginary part of this quantity, 
despite the cr meson has a width larger than 200 MeV for pion 
masses up to around 300 MeV, as shown in Fig. 9, is much 
smaller than its real part, which makes its interpretation eas- 
ier. In the lower panel of the same figure we depict the real 
(blue solid line) and imaginary (red dashed line) parts of the 
quadratic scalar radius of the cr, (r 2 )^. It is notorious that in 
most of this range of pion mass values the square root of (r 2 )^ 
is around 0.5 fm quite independently of the width of the cr 
meson. 

For the physical pion mass we find the values 

(fY s = (0.19 + 0.02) - /(0.06 ± 0.02) fm 2 , 
VV>J = (0.44 + 0.03) - i (0.07 ± 0.03) fm , (93) 
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with the errors calculated as explained in Sec. III. This value 
is almost the same as the corresponding quadratic scalar ra- 
dius for Kit, (r 2 )^ = 0.1806 + 0.0049 fm 2 [126], for which 
the scalar resonance k (or A" ( *(800)), tightly related to the 
cr resonance by SU(3) symmetry [17-21], plays a leading 
role [27, 127]. For comparison, the quadratic scalar radius 
of the pion is {r 2 )" s = 0.65 + 0.05 fm 2 [109]. 17 It is no- 
torious that the value determined for the scalar radius of 
the cr resonance is smaller than that for the pion. It is even 
smaller than the measured quadratic electromagnetic pion ra- 
dius, {r 2 ) n y = 0.439 + 0.008 fm 2 [128]. However, (r 2 )^ is 
similar to the measured K ± quadratic charge radius [129], 
(r 2 )^* = 0.28 + 0.07 fm 2 . Scalar glueballs are expected to 
have even smaller sizes, 0.1-0.2 fm [130, 131]. 

The value obtained for {r 2 )^ in Eq. (93) implies that the 
two scalar isoscalar pions generating the cr resonance are so 
tightly packed that the cr resonance becomes a compact state. 
Whether the two pairs of color singlet valence quarks qq in the 
two-pion state recombine giving rise to combinations of other 
possible QCD states as e.g. q 2 q 2 [17, 41, 132], glueball, etc 
is beyond the scope of our study based on hadronic degrees 
of freedom. In this respect the large Nc evolution of the cr 
pole position [19, 47, 50-52, 133] is enlightening and clearly 
indicates that the cr resonance is not dominantly a glueball 
or a qq resonance. In Refs. [52, 133] it was found that this 
large Nc behavior is compatible with the fact that this reso- 
nance owes its origin to nn interactions becoming a jui reso- 
nance. This large Nc behavior is also compatible with a (qq) 2 
state that fades away as two qq mesons as expected in the 
large Nc limit [84]. This picture on the dynamical generation 
from jtjt interactions of the cr meson is also supported by the 
nontrivial simultaneous fulfillment [133] of semi-local dual- 
ity [133, 134] and scalar, pseudoscalar spectral function sum 
rules [133], both for Nc = 3 and varying Nc- 

On the other hand, for larger values of M„, the cr meson 
closely follows the 2n threshold, as demonstrated in the pre- 
vious section, and its size is then large. Thus, in this range 
of pion masses, the cr meson progressively becomes a two- 
pion molecule and its nature is then much more clear and 
simple (for M„ > 400 MeV it follows from Fig. 19 that 
s/(r 2 ) > 1.5 fm). 18 This can also be related from the behavior 
of the quantity g 2 dG/ds evaluated at s = v (and G evalu- 
ated in the Riemann sheet in which the pole appears). This 
quantity is close to one for a composite meson [ 141-146]. We 
have checked that for the large values of M n in which the cr 
meson is a bound state, we have g 2 dG/ds > 0.8, which points 
to a molecular nature. For values of the pion mass close to the 
physical one, we have instead g 2 dG/ds 0. 

Another interesting point is to consider the relation between 
Fa-(0) and the derivative of the cr pole with respect to the quark 



A recent lattice QCD determination [114] gives (r l )" s = 0.617 ± 0.079 ± 
0.066 fm 2 , or, adding the errors in quadrature, (r 2 )* = 0.6 ±0.1 fm 2 , in 
good agreement with the value given in Ref. [109]. 

A similar value was obtained for the size of the A(1405) resonance in 
Ref. [135], which is also a resonance that qualifies as dynamically gen- 
erated form the meson-baryon interactions [24, 136-139]. In Ref. [140] 
the matter or scalar form factor for this resonance was studied. 
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FIG. 19. Top: the square root of the quadratic scalar radius of the 
cr as a function of M„ is shown for < M n < 600 MeV. Bottom: 
the quadratic scalar radius is represented in the range < M n < 
350 MeV. In both panels, the (blue) solid lines represent the real 
part of each quantity, whereas the (red) dashed line is the imaginary 
part. The points over the curves represent our results for the physical 
case with their statistical errors, Eq. (93). Due to the scale used they 
cannot be appreciated in the upper panel. 



mass. According to the Feynman-Hellmann theorem [147], 
one has the relation: 



dM 2 



FAQ) 
2B 



(94) 



Notice that /v(0) is proportional to B and precisely their ratio 
is not ambiguous. On the other hand, dsa-ldM 2 is given in 
Eq. (31). Then we can write: 



gl(M 2 ) 



dSa- 

dM 2 ' V(V,M 2 ) 
dM 2 

X — 

dM 2 



^ (V( v, M 2 ) - V(s a , MlfGni V, M 2 )) 

(95) 



The dependence of M 2 on M 2 is worked out up to O(M^) in 
Eq. (11) from where one obtains: 

dM 2 . M 2 



dM 2 



1 



(96) 



We show our results for -F lT (0)/2B at NLO and compare 
them with dsa-fdM 2 in Fig. 20, so as to check Eq. (94). In 
the upper two panels we show the real part and in the bottom 
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FIG. 20. Feynman-Hellmann theorem: comparison between 
dsu/dM 1 and -iv(0)/2B, Eq. (95), as a function of the pion 
mass. The (blue) thick solid lines correspond to ds iT jdM 2 at 
NLO, Eq. (31), whereas the (red) dot-dashed lines are evaluated 
at LO, where ds^/dM 2 = dsa-ldM 2 . The (green) dashed lines are 
-iv(0)/2B. From top to bottom, in the first panel the real part of the 
quantities are represented in the range M n , 50 < M n < 600 MeV. In 
the second panel, the same is shown for 50 < M n < 300 MeV. The 
bottom panel shows their imaginary part. 



one the imaginary part. The agreement is certainly remark- 
able for M n < 300 MeV, at the level of just a few percents 
of difference. This range of pion masses is highlighted in the 
second and third panels, from top to bottom. Let us note that 
in Eq. (94) we are comparing two quantities that are obtained 
from the chiral expansion of two different interacting kernels. 
The expansion is not performed on the full amplitudes and 
this is why there is not a perfect agreement, as it is the case 
in the standard perturbative calculations of ChPT [7, 8]. In 



our case the factor V(s lT , Mj) multiplying G//( M„) in the 
right-hand side of Eq. (95) is equal to (V2 + V4) 2 , while O4 
from Eq. (84) only contains Vf, because it is a ChPT one- 
loop calculation at 0(p 4 ). 19 Thus, the differences correspond 
to higher order terms in the calculation of /v(0), beyond the 
0(p 4 ) or NLO calculation of the kernel W(s, s', q 2 ), Eq. (85), 
performed in the present work. 

Another point also worth mentioning is the fact that the left- 
hand side of Eq. (95) does not involve any contribution with 
pion propagators in the external legs but the derivative acts on 
the vertices and intermediate pion propagators in loop func- 
tions. This is also the case in F a (0) once the pion propagators 
in the external legs are removed as explained in Sec. IV C. 



VI. SUMMARY AND CONCLUSIONS 

In this work we have discussed the nature of the cr reso- 
nance (nowadays also called /o(600) in the PDG [1]) by eval- 
uating its quadratic scalar radius, (r 2 )^. This allows one to 
have a quantitative idea of the size of this resonance. 

There are many studies since the nineties based on supple- 
menting Chiral Perturbation Theory with non-perturbative S - 
matrix methods, that clearly indicate a dynamical origin for 
the cr resonance due to the isoscalar scalar 7T7t strong self in- 
teractions [3-6, 12-14]. More recent studies based on the de- 
pendence with Nc of the cr pole [19, 39, 50-52] also corrob- 
orate that this resonance cannot be qualified as a purely qq 
or glueball, with the pole trajectories compatible with the ex- 
pectations for a meson-meson dynamically generated object 
or a four-quark state. In the large Nc limit it is well known 
that loops are suppressed so that the 7T7t rescattering vanishes 
away and then the cr resonance pole disappears according to 
Refs. [19, 52, 133]. These results have been strongly sup- 
ported recently [133] by the simultaneous fulfillment of semi- 
local duality [133, 134] and scalar, pseudoscalar spectral sum 
rules [133], both for Nc = 3 and varying Nc- 

The next question is whether the two pions are loosely dis- 
tributed, so that the cr meson might be qualified as molecu- 
lar or, on the contrary, they overlap each other giving rise to 
a compact object of a size comparable or even smaller than 
that of its constituents. A proper way to answer this ques- 
tion is to determine quantitatively the size of the cr resonance. 
For that we calculate in this work the quadratic scalar ra- 
dius of this resonance obtaining the value {r 2 )^ — (0.19 ± 
0.02) - /(0.06 ± 0.02) fm 2 . Despite the cr has a large width 
the resulting value for the quadratic scalar radius is almost a 
real quantity, which makes easier its physical interpretation. 
This value is very close to the Kn quadratic scalar radius, 
(r 2 )^ = 0.1806 ± 0.0049 fm 2 [126], similar to the measured 
A"* quadratic charge radius [129], {r 2 )f = 0.28 + 0.07 fm 2 , 
and smaller than the quadratic scalar radius of the pion, 
{r 2 ) n s = 0.65 ± 0.05 fm 2 [109]. This means that the cr is cer- 



9 Notice that the derivative with respect to M 2 of the function G(s,M 2 ) is 
proportional to Cq(s, s,M 2 ). 
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tainly a compact object. The square root of its quadratic scalar 
radius is vV>'J = (0.44 ± 0.03) - i (0.07 + 0.03) fm. 

We have further tested our result for the size of the cr by 
considering the dependence of (r 2 )^ on the pion mass. As M„ 
rises the <x meson mass follows the 2ji threshold. This fact 
has been recently observed in the lattice QCD calculation of 
Ref. [124], and was pointed out much earlier in Refs. [19, 20] 
as well as in the more recent work [122]. In such situation, 
with a small binding energy, the expected size of the cr res- 
onance should be definitely larger than that of a hadron. We 
obtain a quadratic scalar radius that increases rapidly as soon 
as the width of the cr meson tends to vanish, which for our 
present NLO fit occurs for pion masses above 330 MeV. In 
this way, already for pion masses around 370 MeV, -<J(r 2 )f 
is larger than 1 fm and diverges for M„ 470 MeV, pre- 
cisely the value at which the cr resonance becomes a zero 
binding energy bound state. In this case, a molecular or nn 
bound state image is appropriate for the cr meson. For even 
higher pion masses, the binding energy is still small which 
gives rise to large sizes for the cr. Nevertheless, we observe a 
steady (albeit weak) tendency to increase the binding energy 
for higher pion masses so that its size tends to dismiss pro- 
gressively, but for the mass range explored in this work it is 
always > 1 .5 fm. The clear tendency of the cr resonance to 
follow the two-pion threshold is a manifest indication for this 
resonance being a meson-meson dynamically generated one. 
For smaller pion masses between 50 and 300 MeV the square 
root of the quadratic scalar radius of the cr meson is rather 
stable with a value around 0.5 fm, independently of its width. 

The value of the scalar form factor of the cr resonance at 
q 2 = 0, Fo-(0), is related via the Feynman-Hellmann theorem 
with the derivative of the cr pole position with respect to the 
pion mass. Within uncertainties, we have checked the fulfill- 
ment of such relation. 

We have studied nn scattering in NLO S U (2) Unitary Chi- 
ral Perturbation Theory as well. We obtain a good reproduc- 
tion of nn phase shifts for 7 = and 7 = 2, and also for lattice 
QCD results of the 7 = 2 scattering length ai and F„. We have 
offered a detailed comparison between different precise deter- 
minations in the literature, including our present calculation, 
of the cr meson mass and width, and of the threshold parame- 
ters a®, b° Q . The resulting average values are a9 = 0.220±0.003 
and b° Q M 2 = 0.279 ± 0.003. For the cr meson pole parame- 
ters we take the mean of the different values with the result 
Ma- = 458 ± 14 MeV and T^/2 = 261 + 17 MeV. Our own de- 
terminations obtained here at NLO in Unitary ChPT are cP Q 
0.219 + 0.005, b° Q M 2 = 0.281 ± 0.006, M lT = 440 + 10 MeV 
and T^/2 = 238 ± 10 MeV. 

ACKNOWLEDGMENTS 

JAO would like to thank Eulogio Oset for an interesting 
discussion. We also are pleased to Guillermo Rios for having 
provided us with the IAM results of Fig. 9. This work is par- 
tially funded by the grants MICINN FPA201 0-1 7806 and the 
Fundacion Seneca 11871/PI/09. We also thank the financial 
support from the BMBF grant 06BN41 1, the EU-Research In- 



k k 




FIG. 21. Diagrams for the one-loop functions Aq, Bq and Co (from 
left to right and top to bottom, respectively). 



frastructure Integrating Activity "Study of Strongly Interact- 
ing Matter" (HadronPhysics2, grant n. 22743 1) under the Sev- 
enth Framework Program of EU and the Consolider-Ingenio 
2010 Programme CPAN (CSD2007-00042). 



Appendix A: Loop functions 

In this Appendix, we give the loop functions used through 
the paper. We start with the scalar one-, two- and three-point 
one-loop integrals, denoted by Aq, Bq and Co, respectively, 
and depicted in Fig. 21. The vector and tensor integrals, de- 
fined later, can be cast in terms of the former. Special attention 
is dedicated to the case of the three-point function, whose cuts 
are also calculated since they are needed in order to evaluate 
the scalar form factor of the cr meson. Notice that all the in- 
ternal masses are equal, as we only have pions as degrees of 
freedom. For this reason, we do not include the dependence 
on the internal mass M 2 in the following (except for the case 
of the function Ao, which does not depend on any external 
momenta). 

Scalar loop integrals 

The simplest one is the one-point loop integral, given by: 

, Cd 4 k 1 
A (M 2 ) = if - . . (Al) 

In dimensional regularization, it results: 

MMl) = tB\ r + log 7?) ' (A2) 

with 

R = ^ 4 - (1 +T'(1) + lo g 47r)J , (A3) 
and « — > 4. 
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The two-point function is: 



BdP 2 ) 



J(2n 



) (k 2 - M 2 + ie){{k - P) 2 - M 2 + ie) 



(A4) 



and analogously it is evaluated in dimensional regularization, 
with the result 



B (P 2 ) 



1 



I6n 2 



M 2 , 

R + log — - 1 -cr(P 2 )log 



cr(P 2 ) - 1 
cr(P 2 ) + 1 



(A5) 

with cr(f 2 ) = yT -AM 2 IP 2 . Since the function is divergent 
and yu-dependent, we define the subtracted function, Bq(P 2 ), 



M 2 

R + log — 

B (P 2 ) = B (P 2 ) - 



I6n 2 



(A6) 



This is the function that will appear in the amplitudes, be- 
cause the term subtracted cancels out with the alike terms in 
the loops and the chiral counterterms. The same procedure, 
applied to Aq(M 2 ), gives 



The usefulness of Eq. (A9) lies in the fact that it is well suited 
for its analytical continuation to the complex plane, which is 
needed in our case, since the cases p 2 — p\ — s a are studied. 
Notice that the residues of the integrals when z — > Z; are zero 
because of the form of the numerators, and also that, since 
z(z - 1) < for z € [0,1], the arguments of the logarithms do 
not cross any cut. 

On the other hand, since the pole of the <x resonance appears 
in the unphysical Riemann sheet, we need to calculate the am- 
plitudes ttttH — > nn in this sheet. This involves the function 
Co(s, s, q 2 ) in this sheet, 20 and this is not so trivial as in the 
case of the function G(s) (see Eqs. (23) and (24)). For that 
purpose, we calculate the discontinuity along the unitarity cut 
of the function Co(s, s, q 2 ): 

ACq = Cq(s + ie, s + ie, q 2 ) - Cq(s — ie, s - ie, q 2 ) , (A15) 

for s > 4M 2 . This can be obtained directly from the integral 
representation in Eq. (A9), and the result depends on the value 
of q 2 . We are interested mainly in the case q 2 < 0, and we find 
two cases: 



A Q = A (M Z ) - 



M 1 



M 1 



167r -(* + log-| = 0, (A7) 



this is why in the amplitudes of Appendix B there is no de- 
pendence on Aq. 

The three-point function is defined by: 

C ( Pl ,p 2 ,q) = iJ — 



A a C = 



A b C = 



1 



r) 4 k 2 - M 2 + ie 



(A8) 



4nA(s 


s, q 2 y 2 

i 


4nA(s 


s,q 2 ) l l 2 
i 


8ttA(s, 


s, q 2 yi 2 



Z- - Z\ 



for q 2 < q 2 an , 



(A16) 



log 



1 - Zl Z- - Zl 



-Zl Z+ - Zl 
1-Z3 



1 



for q 2 an < q 2 < 4M 2 , 
(A17) 



{(k - pi) 2 - M 2 + ie)((k - pi) 2 - M 2 + ie) where we have defined z± = \{\ + cr(s)). In the previous 



and it depends on the three scalars p 2 , p\ and q 2 = {p\ — P2) 2 . 
It is finite and after some manipulations [148], it can be cast 
in the integral form: 



Co(p\,p\,q 2 ) 



1 



1/2 



\6n 2 A{p 2 v p\,q 2 ) 

log f{p\,z) -log f{p\,zi) 
dz + 



X' 

Jo 



dz 



dz 



where we have defined: 



P 2 z(z 



z 


-Zl 




log f{p\,z) 


-loi 


'J{p\,Z2) 


z 


-Zl 




log f(q 2 ,z) 


- log 


f(q 2 ,Z3) 


z 


-Z3 




-l) + M 2 - 


ie , 





equation q 2 n is the so called anomalous threshold, given by 
M 2 q 2 n - -s(s - 4M 2 ) < 0, where the last inequality fol- 
lows from s > 4M 2 . The case that connects continuously with 
q 2 = corresponds to A&Co, which is the required one in 
the calculation of the quadratic scalar radius. For more details 
on the analytical extrapolation chosen see the discussion in 
Sec. IV C. 

As a cross-check of the validity of our procedure, let us 
note that the function Co is related for q 2 — to the deriva- 
tive of the function G(s) with respect to M 2 , denoted by 
dG(s, M 2 )/dM 2 = G. Indeed, one has 



f{p\z) 

A(a, b, c) — a 2 + b 2 + c 2 — 2ab — 2bc - 2ac , 
1 



1 + 



^2 = ~ 



pj-pl-q 2 " 
A{p\,p\,q 2 yi 2 
-2 \ 



1 + 



2 2 
Pi' Pi 



q 



Z3 



1 + 



A(p\,p 2 2 ,q 2 yi 2 ) 

q 2 -p\-p 2 i ) 

A{p 2 ,p 2 ,q 2 yi 2 1 



(A9) 

(A10) 
(All) 

(A12) 
(A13) 
(A14) 



2C (s, s, 0) = G(s) = 



1 



log (A18) 



8n 2 so-(s) b cr(s) + 1 



In the calculation of the scalar form factor Co appears, while 
in the derivative of the cr pole position, v, one has G. Both 
are related through the Feynman-Hellmann theorem, Eq. (94), 
and thus, the unphysical Riemann sheet for the function Co 
must be related to that of the function G(s). When the pole is 
in the unphysical Riemann sheet, we have: 



G„(s) = G,(j) - AG(a) = G,(s) 



4nscr(s) 



(A19) 



20 Notice that we have already taken that s = s', since this will be the case in 
the cr form factor, s = s' = V- 
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If we now calculate the value A/ 7 Co(s, s, 0), we find: 

i 



AfcC (s, s,0) = 



8nscr(s) 



(A20) 



so that 



Co ; //(s, s, 0) = C 0; /O, s, 0) - A 6 Co(*, s, 0) 

i 



= Co ; /(s, s, 0) - 



(A21) 



87ricr(i) 

which implies, as stated, 2Co ; //(s, s, 0) = Gn(s). 

Vector and tensor loop integrals 

Vector and tensor loop integrals appear throughout the am- 
plitudes in Appendix B. We reduce them to the scalar ones by 
means of the Passarino-Veltman method [149]. We start with 
the two-point vector and tensor integrals, defined by 



e)((k - P) 2 - M 2 + ie) 



(A22) 

On Lorentz-invariance grounds, we can write 

Bf = -P"B n (P) , (A23) 

where the minus sign is introduced for convenience, and, per- 
forming the contraction P^B^,, it can be shown that: 

1 



(A24) 

Analogously, the tensor integral can be decomposed as: 

B MV = P»P Y B 20 (P 2 ) + gllv P 2 B 2l (P 2 ) , (A25) 

and the following results, by the appropriate contractions, are 
obtained: 



P 2 B 20 (P 2 ) = 



A (M 2 ) M 2 -P 2 2 1 



(A26) 

, , A (M 2 ) 4M 2 -P 2 , 1 / , P 2 
piB^) = -±-2 + -^B^ - _ \M 2 - — 

(A27) 

For the three-point vector and tensor integrals, we define: 



M ~ J(2nfk 2 -M 2 + k 



(A28) 



[(k-pi) 2 -M 2 + ie)((k - p 2 ) 2 -M 2 + ie) 



and 



C = -tfCn - fiC l2 , (A29) 

CT = p»p\C n + pfrZCa + (fol + plrt) C 23 + , 

(A30) 



where for simplifying the writing we have omitted the argu- 
ments in Cij{p 2 , p 2 v q 2 ). The results for these functions are: 

Cii = (plB (pl) - pip 2 B Q (p 2 ) - (pj - pip 2 )B Q {q 2 ) 

-pl(p 2 l -p l P2)C (plpl,q 2 ))/(2det>H) , (A31) 

C12 = [p\B {p 2 ) - pip 2 B (pl) - (p\ - p\p 2 )B (q 2 ) 

-p 2 l (p 2 2 -PiP2)C (p 2 l ,piq 2 ))/(2det>H) , (A32) 

1 M 2 2 2 2 

+ l - (B (q 2 ) + p\d i + p\C l2 ) , ( A33) 

C 21 = (p\R a - PipiRc) I defH , (A34) 

C 22 = (p\R d - PiPiRb) / dettt , (A35) 

C 23 = (p\R c + P 2 2 Ri, - PiPz{R a + Rd)) 1(2 det-K) , (A36) 
with 



P\ P\P2 ^ 
PlP2 p\ ) 



and 



Ra = 


\B Q {q 2 ) - 


1 V 


~ C 2 4 , 


R b = 


\B (q 2 ) - 


\p\Cn 


- \b (pI) , 


R c = 


\B Q {q 2 ) - 


1 2 r 


- \b*(p\), 


R d = 


\B (q 2 ) - 


■jP 2 C\2 


- C 2 4 . 



(A37) 

(A38) 
(A39) 
(A40) 
(A41) 



Analogously to the scalar loop integrals, we define the sub- 
tracted functions By and Cy by substituting in their expres- 
sions given above Aq — » Aq and Z?o — * Bo- The amplitudes 
nnH — > nn in Appendix B are then written in terms of finite 
and scale independent functions. 



Appendix B: nn H — > nn amplitudes 

In this Appendix, the amplitudes im H — > nn are given for 
completeness. We follow the nomenclature given in Fig. 15. 
The finite contributions to every amplitude are given once the 
infinite and scale dependent terms are cancelled among them. 
In this way, the amplitudes are written in terms of the finite 
and scale independent constants l, as well as the subtracted 
loop functions defined in Appendix A, Bo, etc... The diagrams 
denoted in Fig. 15 by VI, XI and XII, both in the case of 
7r 7r° H — > 7r°7r° and 7r°7r° H — > n + n~, are proportional to the 
tadpole function, Aq(M 2 ), so they do not contribute to the fi- 
nite amplitude, as explained before. 

In the subsequent, unless the opposite is stated, the sub- 
script ; = 1 , . . . , 4 indicates the pion leg with four-momentum 
Pi to which the scalar source is attached. The functions D,, 
corresponding to the inverse of the pion propagators when 
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the scalar source is attached to the ; t h external pion leg, are 
used through this Appendix. These functions were defined in 
Eq. (61). 



b. 7r°n" H -> n + n- 



1. Diagrams I 



a. n°n°H^7i°n a 



6B 

n 



3BM 2 
FiAn^ 



b. ttV H n + n~ 



j^(NLO) 



2B 

B M 2 



(Bl) 
(B2) 



(B3) 
(B4) 



The NLO result corresponds to the LO diagram I multiplied 
by 25Z, with the latter given in Eq. (11). In addition M 2 , F 2 
are expressed in terms of the physical values according to the 
expansions of Eqs. (11) and (12). 



2. Diagrams II 

a. aWtf-mV 



m 



, Tn , 2B s 

y(LO) _ 

1 ~n 


-M 2 

£>i ' 




nrY . 2B s 
2 " F 2 


-M 2 

D 2 ' 




, Tn , 2B S* 
3 ~ Fl 


-M* 
O3 ' 




, Tn , 2B s' 
4 -pa 


-Ml 
D A 




B M°- 4/4(5 


-mD- 


hMl 




16tt 2 




B M 2 4/4(5 


-mD- 


hMl 


FlD 2 


16n 2 




B M 2 n 4/4(5' 


-mD- 


hMl 


^3 


16n 2 




B M 2 4/4(5' 


-mD- 


hMl 


*1*>4 


16n 2 





(B7) 



(B8) 



The NLO results are obtained by multiplying the LO ones by 
36Z and with M 2 , F 2 re-expressed in terms of the physical Ml 
and F%, respectively, according to Eqs. (11) and (12). Notice 
that in addition to the factor Z 2 from the wave function renor- 
malization of the external pion legs there is an extra factor Z 
from the renormalized pion propagator, Eq. ( 10). 



3. Diagrams III 



1 3 



1 A 



j.(LO) _ 



rr 1 = 1 + 



y(LO) 
1 3 



y(LO) 
1 4 



IB 

n 

IB 

n 

IB 

n 

IB 

n 



= =: 1 + 



1 2 



B 

n 

B 

n 

B 

F4 

B 



Ml 
D 2 
Ml 

B>3 

D A 



UMl Ml \U - 3h 



4n 2 
UMl 
4n 2 
l 4 Ml 



Di I6n 2 
Ml 4/4 - 3/3 
757 16tt 2 

4/4 - 3/3 



4tt 2 D 3 16tt 2 
UMl | 4/4 - 3/3 
4tt 2 + £> 4 16tt 2 



The diagrams III and higher in numeration are purely NLO 
contributions. To simplify the writing we then omit the super- 
script NLO in the corresponding amplitudes. 



(B5) 



(B6) 



3B / 4 2 



b. TtV H -> TT + JT- 



(B9) 



(BIO) 
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Ti 

T 2 
T 3 
T 4 



4. Diagrams IV 

a. n n H^n°n 



B l 4 q 2 - hM 2 „ 
Ft 

B J 4 q 2 - l 3 M 2 n 
Ft 

B J 4 q 2 - l 3 M 2 n 
Ft 

B Uq 2 - l 3 M 2 n 
Ft 8^2 





M 2 n 


1 + 






D x 




Ml 


1 + 






D 2 




Ml 


1 + 




D 3 




Ml 


1 + 









Ti=- 
T 2 

T 4 



B 



24n 2 Ft 
B 

'24n 2 Ft 
B 

TAn^Fi 
B 



24n 2 Ft 



hs + W + u) - 2M„(h + 2l 2 ) - 
h s + l 2 {t + u ) - 2M 2 n (h + 2/ 2 ) - 
j ( h s' + ~W + u ) - 2Ml(h + 2l 2 ) - 
h s' + T 2 (t + u) - 2M 2 n (h + 2l 2 ) - 



r(s , t , 


U ) 


£>i 




P(s , t , 


u ) 


D 2 




P{s',t', 


u ) 


D 3 




P(s',t , 


u') 


D A 



(B15) 



(Bll) 



Ti = 
T 2 
T 3 
T 4 



b. ttV H 



B kq 2 -l 3 Ml Is -M 2 n 



Ft 8tt 2 \ D\ 
B ~kq 2 -hMl Is - Ml 



Ft 8tt 2 \ D 
B hq 2 -hMl Is' -Ml 



Ft 8tt 2 \ D 3 
B hq 2 -hMl Is' -Ml 



Ft 8n 2 



D 4 



5. Diagrams V 

a. ;rVtf->»rV 



(B12) 



6. Diagrams VII 

a. Wfl^A 



B ( , , 2o 2 - Mi \ - 

T\ = -jt \lq 2 + D l+ M 2 K H - " \ B (q 2 ) 

B f , 9 2o 2 -M 2 ; 
T 2 = - — \2q 2 + D 2 + M 2 n H D2 ' \ , 

B f , , 2o 2 - Mi \ - 

T 3 = -j i ^2q 2 + D 3+ M 2 n H ^ " \ B (q 2 ) 

Bio 7 2q 2 -M 2 \ 

TA = ~y A W + D A + Ml q ^ " \ B (q 2 ) 



(B16) 



r, = 



T, =- 



T t 

Ti 



B 


h + 2l 2 


Ft 


24tt 2 


B 


h + 2/ 2 


Ft 


24tt 2 


B 


h + 2/ 2 


Ft 


24tt 2 


B 


h + 2l 2 


n 


24-7T 2 



D\ + 2Ml - 



D 4 + 2M; ■ 



s 2 


+ t' 2 


+ u' 2 


-4Mf 






Di 




S 


+ f 2 


2 

+ u 


-4Ml 






D 2 




s' 2 


+ t' 2 


+ u 2 


-4Mt 






D 3 




s' 2 


+ t 2 


+ u' 2 


-4Mt 


Da 



(B13) 



b. H 7i + n- 



In these amplitudes, we define: 

P(s, t, u) = 4M% + 2l 2 ) + hs 2 + hit 2 + u 2 ) 

-2M 2 K {2hs + 2i 2 (t + u)) (B14) 



Ti = 
T 2 = 
T 3 = 
T 4 = 



b. n Q n° H ^ tt + 7t~ 



B , n f 2q 2 - Ml ) - , 

(*-M*){l + ^}so(^ 

(s'-Ml){i + - 2 ^-)B 0iq 2 y 



B 

n 

B 

n 

B i , 2 \{ 2q 2 -M 2 \ . 

(s'-Ml)\l+ - * \B (q 2 ) 



Di 



Ft 



Da 



(B17) 
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7. Diagrams VIII 

a. n n H^n°n 



^-channel diagrams: 



Ft 2 



T 2 = 
T 3 = 

r 4 = 



Ft 2 



Ft 2 



n i 2 

f-crossed diagrams: 

2B [Ml 
Fx =-— <!^ L + 

r 2 = 
r 3 = 
r 4 = 



2 



2 



F* 2 



Fl| 2 



w-crossed diagrams: 



Ft 2 



To = - — r 



2 



r 3 

T 4 





2s Ml 


Di 




\M\ + s 2 - 


2s Ml 


u 2 




\M A n + s' 2 - 


2s'M 2 n 


n, 




\M\ + s' 2 - 


2s'M 2 n 


U A 




\M A n + 1' 2 - 


2f'M 2 ' 


Di 






2f M 2 






\M\ + 1'~ - 


2f'M 2 






3 n *A .9 

+ f - 


2t Ml 


n. 




|M* + m' 2 - 


2u'M 2 n 






|M4 + « 2 - 


2u Ml 






\Mt + u 2 - 


2u Ml 



B (s ) , 

BoO ) , 
B (s') , 

B (s'). (B18) 



B (t') ■ 
B (t) 
Bo(t') . 
B (t) 



2B(Ml \M A n + u' 2 -2u'M 2 n 

"n\~T + d~ 4 



B (u') , 
Bq(u ) , 
Bq(u ) , 
B (u') . 



b. ttV H -> n + re- 
in the t- and M-channel amplitudes, we define: 



(s + u- 2M 2 )(M 2 n -t)+ M 2 t - 



(B19) 



^-channel diagrams 
B 

T\ = -(s — i 

71 

To — t(s — j 

n 

B , 

n 

1 IT 

f-crossed diagrams 



s 


+ Ml 




Di 


s 


+ Ml 




F>2 


s' 


+ Ml 




D 3 


s' 


+ Ml 



D 4 



(B22) 



T x = 


2B Q(s ,f. 


,«') 


n Di 




T 2 = 


2B Q(s,t, 


,u ) 


n d 2 




T 3 = 


2B Q(s',t', 


u) 


n °i 




Ta = 


2B Q(s',t . 


,«') 


n °4 





(B23) 



M-crossed diagrams: 



Tx = 


2B Q(s ,u' 


\o 






T 2 = 


2B Q(s ,u 


,t) 


n D 2 




T 3 = 


2B Q(s',u 


,f) 


n ° 3 




Ta = 


2B Q(s',u' 


\t) 


n d a 





(B24) 



8. Diagrams IX 

In these amplitudes, the scalar source s can be attached to 
one of the two four pion vertex (recall Fig. 15), which we 
denote here by the subscript ; = 1,2. 



a. n° n° H -> ;rV 



(B20) i-channel diagrams: 



T x = (2s +Ml)Bo(s) 



B 



T 2 = —(2s'+Ml)B (s'). (B25) 



f-crossed diagrams: 



T x = —,(2t +Ml)B {t) 



n 

B 



+ t(s + u- 4M 2 )B 2 o(0 + 2t(u - 2M-)B 2x (t) (B21) 



r 2 = — (2t'+Ml)B (t') 



(B26) 
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w-crossed diagrams: 



T-i = j- A {2u + M 2 )b ("), 



b. ttV H tt+tT 



^-channel diagrams: 



Ty =-^(4* -3il^)5 (j) 



r 2 = -^(4/-3M^)5 (*') 



r-crossed diagrams: 



7! =-— (r -2M 2 )fi (O 



T 2 = -—(lf -2Ml)B (f). 



M-crossed diagrams: 



Ti 
T 7 



= ~(« -2^)B (") 



9. Diagrams X 



^-channel diagrams: 



T = ~ | (2(5 - M 2 .)^' - M 2 ) + A<) C (J, q 2 )+ 
+ Ml (5 (s) + BoO') + 2B (^ 2 )) + * + " ~ q B (q 2 ) 



w-crossed diagrams: 

T = ~ I (2(« - M 2 )(m' - M 2 ) + M*) C (u, u', q 2 )+ 



T 2 = ±. (2 U ' + M 2 ) fl o («0 ■ (B27) + M * ( 5 o(«) + *<)("') + ^ (q 2 )) + 



2 \ U + u' — q l 



Mr) 



(B28) 



(B29) 



(m' - 2M 2 ) B (u') . (B30) 



(B31) 



f-crossed diagrams: 



T = ~ | (2(f - M 2 )(f' - M 2 ) + M*) C (t, f, q 2 )+ 



i 2 

M 2 (5 (0 + SoCO + 2Bo(? 2 )) + /+ ~ g 



(B32) 



(B33) 



s-channel diagrams: 
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r = -^4 IO' " <)Co(^, s',? 2 ) + (* + ✓- 2Ml)B {q 2 )+ 



+ (s' - M 2 n )B (s') + (s- Ml)B Q ( S ) 
r-crossed diagrams: 



(B34) 



2fi/(2M 2 -f)(2M 2 -f') 

2Mz - 1 , - , _ x 
— ^— (B (q 2 ) + B (tj) 



C Q (t,t',q 2 )+ 
2Ml - f 



, t + t'-q 2 - 2n ( s - s') 2 - (u - u') 2 - 
+ ^ #o(<7 ) + r 

- (i + 5 1 ' — U — U 1 ) C24(t, t' , q 2 ) 

w-crossed diagrams: 



(B (q 2 ) + B (t'j) 
C23(f,f',^ 2 ) 

(B35) 



^ = 7 



2fi /(2M 2 - m)(2M 2 - «') 



Cq(u, u ', q 2 )+ 



2Ml - u , - , _ x 2M 2 - m' , _ , _ , x 

— | (Bo(<? 2 ) + flo(«)) + — ^ (Bo(? 2 ) + flo(«')) 



U + «' - <? 2 6 2 (a - s') 2 - (t ~ t') 2 
- 4 B (q ) + 

(s + s' — t - t') C2a(u, u, q 2 ) 



10. Diagrams XIII 



a. nVff-^W 



Cn(u, u',q 2 ) 
(B36) 



T= — (25M 2 + 8<? 2 )B (? 2 ). (B37) 



r = — (5(5 + s' - M 2 ) + 9 2 ) B (? 2 ) - (B38) 
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